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CHAPTER 0 : INTRODUCTION 


Magnetic resonance imaging (MRI) is regarded as today’s most powerful non- 
invasive diagnostic tool in medicine. Many areas of science, engineering and 
technology have contributed to its growth and are continuing to do so. This thesis is 
an offshoot of the collaboration between the Indian Institute of Technology, Kanpur 
and the Sanjay Gandhi Post Graduate Institute of Medical Sciences, Lucknow, started 
by the enthusiastic efforts of Professor R.K.S. Rathore, Department of Mathematics, 
I.I.T. Kanpur and Professor R.K. Gupta, Department of Radiology, SGPGIMS. 

One of the achievements of the collaboration has been the observation of 
Professor Rathore that the MRI images look blurred to him, and to de-blur which he 
proposed an effective iterative technique currently being used at the SGPGIMS. The 
technique has been demonstrated at several radiological meetings, including the 5* 
National Symposium on Magnetic Resonance, IIP Dehradun, a brain storming session 
at SGPGIMS, and a National Workshop on MRI Imaging with Emphasis on Neuro- 
Sciences at SGPGIMS, attended by leading MRI experts of the country and abroad. 

The thesis, in brief, is a study of this technique from a classical point of view, 
in which the models of the blurring explored are taken to be the celebrated Fej6r sums 
<Tn, the generalized Jackson operators Ln^, and the variation diminishing integrals Vn 
of De la Vall6c-Poussin. The operators < 7 „, L„,^, and V„ are linear positive operators, 
well known in the theory of approximation of functions. Of these, the last two 
produce blurrings of the respective orders tWjf/d/n) and which are 

precise. However, for the Fej6r sums a„, no such general comprehensive order is yet 
known (they blur Lip a functions with order l/n®, 0 < a< 1, and Lip 1 functions with 
order (In 72)/n). These three distinct behaviors have been the chief reason motivating 
their inclusion in this study. In the thesis we obtain theoretical direct and inverse 
results about the related sequences of de-blurrings and also consider their application 
on actual MRI images, with the numerical simulations showing their effectiveness. 

This introductoiy chapter describes some of the basics of MRI, the de-blurring 
technique of R.K.S. Rathore, an overview of the main theme of approximation theory 
that degree of approximation depends on the smoothness of the approximated 
function, followed by a summary of contents of the rest of the chapters of the thesis. 



0.1 BASICS OF MR! 


For the basics of Nuclear Magnetic Resonance (NMR) imaging, a number of 
excellent texts, e.g., T. Farrar and E. Becker (1971), P.T. Callaghan (1991), P. 
Mansfield and P. Morris (1982), C.L. Partain, R.R. Price, J.A. Patton, M.V. Kulkami 
and E.V. James Jr. (1988) are available. [References are arranged alphabetically by 
authors names (year of publication) at the end of the thesis] 

A schematic diagram of an NMR imaging system is shown in Fig. 0.1. 



Fig. 0.1 ; Schematic diagram of a general NMR imaging system 

The patient lies in the probe area. The pulse sequence chosen determines the 
action of gradient controller, transmitter and the receiver coils. The receiver gets 
signals in the audible range which are amplified and digitized for data measurement 
and storage purposes. The signal amplification, measurement and its digitization are 
some of the possible sources of system blurring of interest to us. The CPU, which 
computes the reconstructions, may also be instructed to apply additional filters, 
depending upon the noise characteristics of the pathological region under 
consideration and the number of acquisitions of the data for a particular experiment to 
improve upon the signal statistics. The reconstructed images could be displayed, 
stored and put on a transparency. The physics of MRI and related mathematical 
equations follow next. 










BULK MAGNETIZATION VECTOR 


Following R.R. Price, W.H. Stephens and C.L. Partain (1988), when an 
external static magnetic field B© is applied, the alignment of the randomly oriented 
magnetic moments of individual nuclei creates a macroscopic magnetization M 
pointing along the axis of the external field as shown in the Fig. 0.2. 



Fig. 0.2 Alignment of the individual magnetic moment vectors creates a macroscopic magnetization M 
pointing along the axis of the external static magnetic field Bq 


PRECESSIONAL FREQUENCY 

In the presence of an applied magnetic filed Bo, the magnetic momentum 
vectors experience a torque and process about the axis of the magnetic field with a 
rate given by the Larmor’s relationship 

^271 2n ' 

where /is the resonance frequency, O) is the angular frequency and y is gyromagnetic 
ratio. The gyromagnetic ratio is equal to the ratio of the magnitude of the magnetic 
moment of the nucleus to the magnitude of its spin angular momentum. 


RADIO-FREQUENCY MAGNETIC FIELD 

Magnetic resonance absorption can be detected only if transverse 
magnetization is created, since it is the transverse component Mxy of M, that is time 
dependent and thus according to Faraday’s law of induction, can induce a voltage in 
the receiver coil. 



Referring to Fig. 0.3, when the radio frequency magnetic field Bi is turned on, 
it appears perpendicular both to M and to the angular momentum of M, with the result 
being that M will try to process around Bi with the effect that M tips away from its 
original orientation. When this happens M acquires a component that is perpendicular 
to Bo causing it also try to process around Bq. With Bj constantly changing at an 
arbitrary frequency M just w'obbles about it original orientation. However if the 
frequency of Bi field is equal to the Larmor frequency of the field Bo, due to 
resonance the effect of much smaller Bi is cumulative, causing the angle of the tip to 
continue to become larger and larger. If the duration of the Bj field is such that the net 
magnetization is rotated by an angle 90°, it will become transverse. The angle of 
rotation 0, the RF flip angle, is given by 6 = yBiX, where Bi is the amplitude of the RF 
field and x is its duration. Hence, RF angle can be controlled via the amplitude or the 
duration of the Bi field. 


t 



Fig. 03 (a) Radio frequency filed B] rotating synchronously with the processing spins. 

(b) A 90° flip to magnetization M brings the entire magnetization in the xy plane 


FREE INDUCTION DECAY (FID): 

When the rf-field Bi is removed the magnetization is subjected to the effect of 
the static field only and hence precesses about it. Now by placing a conductive coil 
around the sample, with the axis lying on the xy-plane, it will induce an electric signal 
in the coil This is the only measurable quantity in an MR experiment and is known as 
Free Induction Decay (FID). 




THE BLOCH EQUATIONS 


F. Bloch (1946) explained the motion of the macroscopic magnetization in the 
presence of an applied magnetic field in terms of phenomenological differential 
equations. The classical equation of motion (see, eg., A. Abragam (1961), and T.C. 
Fairar and E. Becker (1971)) of a magnetic moment in a magnetic field B states that 
the rate of change of angular momentum of the spinning nucleus, p, depends upon the 
torque, fx xB , exerted on the magnetic moment by the applied field as: 


Using the fact that p.=7p , 


dt 


= \x xB. 


d\i. 



= 7 p.xB . 


For the macroscopic magnetization M = E p,„ above relationship will be 

^ T> 

=yMxB. 

dt 

In general B consists of both the static applied field Bo, and the magnetic vector of the 
radio frequency (rf) field B]. The rf-field can be thought of as a field rotating in the 
jcy-plane at angular frequency co. During the RF pulse, the magnetic field is 

B(r) = (t) cos CO ti - (t)sin(urj + Bk , 

where B = BQ+h, h being the contribution of the gradients. Thus Bloch equations 
become 

= YM^.B-yM^B, (r)sin cot , 
dt 

dM^. 

--YM.Bi coscot-YM^B, 

dt 

and sin cot - costut . 

dt 

With the initial condition (m, (0) (0) M , (0))^ , the solution of the above 

system of equations is given by: 

\ <2 . . 1 

M^{t) = M^{(i)\ ^ + ~ 2 — 77^°® + -==== sin ct sm cot 1 

a +b fa^+o J 



D 


+ M^(0)\-y4= — s'moxcoso)!- cos ct sin cot 

[^Ja^+b^ I 


, ,, ,rvM b . . , 

+ M . ( 0 )^ — T cos cot — ; cos ct cos cot — — sin ct sin cot ^ , 


a^^b 


a^+b^ 




M At) = M (O)-j-^; — 7vsin(wr + -^; — —cosctsincot — p===rSincrcostyr 


a^+b^ 


a^+b^ 


■^a^+b^ 


J 


+ M,(0) 




■yfa^ +&' 


r sin c? sin fi» + cos c? cos car 
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1 


+ M,(Q) - - - " , {a^ +Z?^coscr), 


a^+b^ 


where a = + ca , b = YB^ and -i-b^ . 


THE SPIN-LATTICE AND SPIN-SPIN RELAXATION 


F. Bloch (1946) assumed that spin-lattice and spin-spin relaxation could be 
treated as first order processes with characteristic times Ti, spin-lattice relaxation or 
longitudinal relaxation time, and T 2 , spin-spin relaxation or transverse relaxation 
time, respectively. The transverse components Mx and My decay back to their 
equilibrium value zero, while the longitudinal component Mz returns to its equilibrium 
value Mo- 

In the absence of RF field and considering Ti, T 2 relaxation, Bloch equations become 


dM^ 

dt 

dM^ 

dt 


= yM B 


M. 


-yM^B- 


'2 

M. 


dM. _ M, -Mq 

T, 


dt 



With the initial condition (a/^(0) M ^ (0) Mj(0))^ , the solution to the above system 
of equations is given by: 

M ^{t) = exp(-? / )'{a/ ^ (0) cos + M ^ (0) sin yfir }, 

My{i) = exp(-r / r, ^ (0) cos yBt - (0) sin ySr}, 

and M^(t) = M^+(M^(0)-M^)exp(-tIT,). 

Thus, if Mt (t) denotes the magnitude of the transverse component of the 
magnetic moment vector and Mi (/) denotes the longitudinal component itself after a 
duration of time t of relaxation, and Mj (0), and Mi (0) denote the corresponding 
values at the start, the Ti and 7% relaxation effects are given by 

M^(t) = M(l - ) + M^(0)e"''^' , 

Mj-it) = Mj(0)e~‘^^^ (true decay), & 

Mt (t) = Mt (0)e~“^^ (FID decay), 

where M = Mq, the equilibrium value, and Tj* denotes the transverse decay constant m 
the presence of possible magnetic field inhomogeneities, as compared with Ti, which 
ought to be a purely atomic constant. 

In the case of the standard spin-echo (SE) experiment (see J.P. Jones (1988)), 
the first starting sequence has the event-wise associated longitudinal and transverse 
magnetic moment components given by Table 0.1: 


Event 

Longitudinal Component 

Transverse Component 

Start 

M 

0 

90° 

0 

M 

TE/2 



180° 

-M(\-e~^^^') 


TE/2 





i— measure 

End 


j^e-rBiT,^-TDir; 


Table 0.1:The First Spin-Echo Sequence 






The second spin-echo sequence, which is used for acquiring the Ti, T 2 and 
density weighted SE-images, has the corresponding longitudinal and the transverse 
components given by the following Table 0.2: 


Event 

Longitudinal Component 

Transverse Component 

Start ! 

= 0 

90° 

0 


TE/2 


M(l-2c 

180° 


M (1 - -H )e-^^''^'^ 

TE/2 


M(1 - 


! 

<— measure 

End 

I 



I 


Table 0.2: The Second Spin-Echo Sequence 

TWO AND THREE DIMENSIONAL IMAGING 

The basic idea of a spatial labeling by an encoding process was given by P.C. 
Lauterbur (1973), who also obtained the first MR-images using the projection 
reconstruction. His technique is known as frequency encoding. 

The Fourier magnetic resonance imaging proposed by A. Kumar, D. Welti, 
and Richard R. Ernst (1975) uses both frequency and phase encoding to obtain the 
Fourier coefficients of the tissue parameter weighted distribution of the magnetic 
moment vector M in the region R of concern. Thus, using the standard spin echo 
sequence the data provided by an MRI scanner consists of the mathematical quantity 
given by 

F{t^,t) = , 

R 



where y is the gyromagnetic ratio; Gx, Gy, G, are the respective x, y, and z gradient 
strengths; ty are the y and z phase encoding times and 5 denotes the effective 
intensity of the signal at the point (x, y, z). Different choice of t, ty, and tz enables us 
to measure the 3D-Fourier coefficients of the signal intensity function. 

In the normal 2D-imaging one simply does not use the z-gradient phase 
encoding to obtain the 2D-Fourier coefficients of the particulai' planar slice under 
consideration. However, the slice selection takes place via a z-gradient on the basis of 
which the appropriate rf-frequency range is invoked. 


TISSUE CONTRAST 

In MR-imaging, the contrast between tissues depends upon local values of 
spin density p, relaxation times Tj, T 2 and inter-pulse time intervals - echo time TE, 
repetition time TR, delay time TD, etc. Using the Bloch equations, signal expression 
for spin echo (SE) sequence (90°-TE/2-180°-TE/2-measure-TD), where TD = TR-TE, 
can be written as (see, e.g., G. Sebastiani and P. Barone (1991) & Tables 0.1-2), 


where is Mo is proportional to the proton density, p of water in the tissue. 


The following Table 0.3 shows different contrasts in the image which can be 
achieved by suitably choosing the sequence parameters TE, TR in the case of spin 
echo SE-puIse sequence. 


I Contrast 

I 


Ti 


T2 


Proton 

Density 


TE -Setting 


TR -Setting 


As short as possible (preferably Comparable to T\’s of two tissues 
less than 20 ms) j (towards the shorter T\ is better) 

Comparable to Tfs of the two j Long compared to the Tfs of the two 
tissues (towards the longer T 2 is \ tissues (2000 ms or greater) 
better) ! 1 

As short as possible (preferably ; Long compared to the Ti ’s of the two ! 
less than 20 ms) ■ tissues (2000 ms or greater) ! 


Table OJ: Different contrasts obtainable in Spin Echo Imaging by varying TE, TR 
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Fig, 0.4: Ti (TE = 15 ms, TR = 300 ms), Tj (TE = 80 ms, TR = 2000 ms), and Proton density (TE = 

20 ms, TR = 2000 ms) weighted SE-images 

In each of the main chapters 1-3, these Tu T 2 , and p-weighted images have 
been used as reference images (or actual function f) for simulations. Subsequently 
they are blurred (Aq/) by the filters concerned and the performance of the de-bluning 
(Ai/, A/, A/,. . .etc.) has been documented in the chapters by computing different error 
measures associated with the degree of approximation {ej =f-Ajf,j> 0). The MRI 
data for the thesis has been acquired on a 1.5 Tesla machine at SGPGIMS, Lucknow. 


THREE MAIN TEST CASE IMAGES 


The following three test case images (see Fig. 0.4) are some Ti, T 2 and Proton 
density weighted images of the same cross-section, obtained by spin echo sequences 
with appropnate TE, TR values. 
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0.2 ON THE ITERATIVE DE-BLURRING SCHEME 
OF R.K.S. RATHORE 


The need for filtering in MRI is different from that in digital image processing. 
In the latter it is to compensate for a degradation that the image itself has suffered, 
while in the former it is to reduce the data noise associated with the high frequency 
terms in the K-space. The filtered image, however, looks blurred or unfocussed. Thus 
assuming that the effect of the noise gets diminished, to achieve a finer diagnosis at 
SGPGIMS, Lucknow, Professor R.K.S. Rathore proposed a de-blurring based on 
some good approximation of a localized spatial blurring model suitable for the MRI 
scanner in use. His spatial de-blurring technique consists of a surprisingly effective 
iterative inverse filtering procedure, through which the noise contribution continues to 
remain dissipated, whereas the image structure tries to re-focus and thereby acquire 
approximately the original sharpness so that the overall visual diagnostic resolution of 
the cross section gets sufficiently improved. 

Following R.K.S. Rathore (1994), let the actual cross section be/and the local 
blurring operator be denoted by U. The filtered image could then be modeled as Uf+ 
N, where N is the residual dissipated noise in the image. Note that, as mentioned 
earlier, U and N arise due to a variety of factors including the signal amplification. 
Postulating an approximation A of U, the proposed class of sharpening operators is 
given by 

'"C/^-',/n=l,2,3,4,... 

If ||zi|| = 11/-AII < 1, Arr^m A-\U-A)f, whereas |1A„N1| < [(l-l|Ain/(l-||A||)]l|Nl|. 
Thus the resulting error on the application of Am for a moderate m would have the 
order of approximation of the modeling error and the noise could increase by 
maximum multiplicative bound (1-1|A||)'‘, which would normally be of the order of 
unity. In the reconstructions shown in this section, A was chosen to be an appropriate 
blurring model (based on a variant of the Fej^r sums) quite suitable for the present 
SGPGIMS MRI scanner. Note that in general, the order of approximation would be 
higher for the smoother structures in the cross section [R.K.S. Rathore (1994)]. 
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It may be noted that most of the usual methods e.g., Wiener filtenng or inverse 
filtenng by dividing the Fourier data by the Founer transform of the point spread 
function (psf) do not give good results in the MRI context. This is because the former 
Itself causes a blumng, while the inverse filtenng leads to severe artifacts in the 
image arising due to meagemess of the higher order frequency content of the psf. The 
problem is further complicated because of the erratic behavior of the signal amplifier 
for the meager K-space signal at the higher frequencies, as example in Fig.0.5 shows: 



Third Iterate A4 Inverse Filtered Image 


Fig. 0 . 5 : Top left - scanner image (Aj), top right - A3, bottom left - A4, bottom right - inverse filtered 

In Fig. 0.5, the scanner image Ai is that of a pig muscle (TE = 80, TR = 2750) 
laden with of cysticercal cysts. The inverse filtered image consists of pure artifacts, 
whereas the resolution enhancement in the iterates A 3 and A 4 is visibly obvious 
allowing a clear counting of the number of cysts from a magnified version. 
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His technique was clinically tested on data from 250 patients who had MRI for 
different pathological conditions involving brain and spine. The images obtained with 
different pulse sequences like conventional spin echo, gradient echo and magnetization 
transfer, were subjected to the technique to look for resolution enhancement. The 
enhancement of resolution was independently confirmed by two experienced MR 
radiologists with more than 10 years of exposure to MR image interpretation. The 
resolution was considered to be enhanced if the lesion became more conspicuous on 
applying the technique than without it. Further test of resolution was done on a pig 
muscle full of cysticercal cysts (a section shown in Fig. 0.5). The number of cysts seen 
before using the technique and after using the technique were compared with actual 
counting on ex vivo and were more on post processed images compared to preprocessed 
images and were found equal to the actual number of cysts. 


Fig. 0.6: Iterative De-bluning in SB and GE-images 

Fig. 0.6 has three rows of images, in each of which the first image Ai is the 
unprocessed one. The first two rows subsequently display the first to third iterates A 2 , A 3 , 
A 4 , while in the last row the second image corresponds to the third iterate A 4 . 
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For the fist routine T 2 weighted, SE-image (TE = 80 ms, TR = 2200 ms), the 
mean (n) and standard deviation (a) statistics were computed for some three ROI’s 
and are displayed in the Table 0.4. It is seen that p in different regions is not affected 
much, while the a increases in the sharpened images with the increase in the index m. 
This is expected, as a higher resolution implicitly means more discrimination and 
consequently a higher variability in regions of a cross-section having less smoothness. 



REGION 1 

REGION 2 

REGION 3 

IMAGE 

P 

B 

P 

a 

P 

a 

Ai 

169.1 

B 

124.7 

B 

63.5 

59.1 

A2 

169.1 

44 1 

125.3 

19.8 

63.7 

61.4 

A3 

169.4 

46 5 

125.3 

22.3 

64.1 

61.8 

A4 

169.3 

47 8 

125 3 

23.9 

64.1 

62 3 


Table 0.4: Mean (ju ) and S.D. (cr) in Different Regions 


Row 2 of Fig. 0.5 corresponds to a Ti weighted, gradient echo BOLD-image 
(TE = 46 ms and TR = 93 ms) and equally well shows the progressive sharpness of 
the resolution in the iterates A 2 , A 3 , A 4 . 

Row 3 of Fig. 0.5 corresponding to a standard T 2 weighted SE image (TE = 80 
ms, TR = 2200 ms) shows multiple lesions in the lateral rectus of the right eye of a 
patient. The eccentrically placed scolex is more clearly visible in the post-processed 
image A 4 compared to the preprocessed Ai. Demonstration of the scolex is essential 
for the diagnosis of cysticercal cyst, as it is pathognomonic of this condition. 

In the numerous reconstructions obtained by his technique, no artifacts or any 
lesion degradations were found. The sharpness achieved in all cases had been 
encouraging throughout. A few follow-ups confirmed the diagnoses retrospectively. 

The technique, thus, was found to be a very promising tool to remove the 
blurring in MRI images. The resulting images did not show artefacts or lesion 
degradation and were sharp. In the clinical practice they have led to several successful 
diagnoses, which otherwise would have at best been categorized as borderline cases. 
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0.3 AN OVERVIEW OF APPROXIMATION THEORY 


An MRI reconstruction is basically a Fourier reconstruction, which is related 
with the area of approximation theory dealing with the approximation of 27T-periodic 
functions by trigonometric polynomials. The main emphasis there is on the 
approximation of 2n:-periodic continuous functions, the class of which is denoted by 
Cjn ■ 

The discontinuities in the functions responsible for the well-known Gibbs 
phenomena, etc. need to be treated by a special asymptotic analysis, distributional 
methods, or by regarding them as points of very rough smoothness. For MRI purposes 
it seems useful to regard a cross section as a sum of several functions with graded 
smoothness. Thus, while one analyses a particular region, because of the localization 
of errors in Fourier reconstructions, one has mainly to keep the smoothness of the 
region under consideration. 


If {La?} is a sequence of bounded linear operators on C 2 n , we shall use the 
writing convention LA/{f, x) = LN(fit); x), where “r” represents a dummy variable. It is 
to be understood in the sense of the following working expressions: 




sin 


(2n + l)(t - x) 


sin 


t-x 


-dr = 5„(/(0;x). 


Several definitions become easy with this usage, e.g., the k-th moment of at x 
could be defined as: fik_N{x) = Lniit-xf; x). 

For approximating functions of two variables, which constitute the MRI 
images, we use the derived Cartesian array defined by 

ify x, y) — - s), x), y) = s'), y), x), 

where Lm acts on the first argument “t” and Ln on the second argument “ 5 ” of fit, j). 
For such Cartesian product operators, writing 

||/|| = supsup|/(x,y)|, 

X y 



16 


degree of approximation results for are immediately obtainable from those about 
Lm and L^, for: 

li-m/ - /II = ||i„ (i-,/) - 4 = 0(«)(i) + ?.(>)), 

is equivalent to 

||i»/ - /II = W -4 ~ , ||i„/ -/I = 0(«)(i)), W ^ ~ . 

Similarly, if w{g ;t) is a measure of smoothness (w(^ ;'r)— > 0,T — > 0), for functions g 
of one variable and for functionsy(x, y) of two variables we define 

w(/;r,T) = supw(/(.,y);T) + supw(/(x,.);0, 

> AT 

then w'ithT/^ — > 0, N and 0, Af — » oo , we have 

Kf-f\\ = 0{wif-,x^)\M^^. & \\Lt,f-f\\ = Oiw{f-,x,)lN-^oo 
are together equivalent to 

In view of this discussion, since the MRI filtering is of the above Cartesian 
nature, the degree of approximation results in the thesis have been derived for the one 
dimensional case only. The numerical simulations, however, correspond to 2-D 
arrays, where a convolution in x-direction is followed up by another in the y-direction. 
In the case of convolutions the operations being commutative: LmL„ = LnLm. 

For literature on basic results in approximation theory dealing with the degree 
of approximation, we refer to the excellent texts of N. Wiener (1932), A. Zygmund 
(1959), I.P. Natanson (1964), A.F. Timan (1966), P.P. Korovkin (1960), G.G. Lorentz 
(1966), P.L. Butzer and H. Berens (1967), H.S. Shapiro (1969), R.A. DeVore (1972), 
and Z. Ditzian and V. Totic (1987), to quote a few. The historical material on the 
Fourier series is very nicely dealt with in the classic text of H.S. Carslaw (1930). 

To summarize, Joseph Fourier in the year 1822 asserted the trigonometric 
series expression of an arbitrary function. The proper systematization began with the 
classical approximation theorem of Karl Weierstrass in the year 1885 stating that an 
arbitrary Czn function can be approximated arbitrarily closely by a trigonometric 
polynomial. 
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The degree of approximation by trigonometric polynomials started with the 
celebrated Jackson theorem, which constituted the doctoral dissertation of Dunham 
Jackson (1912) stating that: 

Enif) < AfiXf; m, n = 1, 2, ... ,/e Cjn, 

where A is an absolute constant independent of/ and n, oXfi d) stands for the classical 
modulus of continuity of/ defined by 

oif\ d) = max|;,|<5 |1 /x+/z)-/(j:)||, 5 > 0, 

and Enif) denotes the error in the best approximation of / by trigonometric 
polynomials of order n. Thus, if oj(f, S) = 0(5®), E„(f) = 0(n'®), 0 < a < 1. The 
estimates of the type of Jackson are referred to in the literature as Jackson type 
inequalities, direct theorems, or a priori error estimates. The proof of Jackson was 
constructed with the help of a sequence of linear trigonometric polynomial operators, 
known as Jackson operators. These correspond to the case p = 2 in the sequence of 
generalized Jackson operators {L„j, } considered in Chapter 2 of this thesis. 

Sergi N. Bernstein (1912) came out with the converse of the above result of D. 
Jackson, stating that En(f) = 0(n'®), 0 < a < 1, in turn implies that off, S) = 0(5“). 
Bernstein, however, could only prove that co(f, 5) = 0(5 In (1/5)), in the case a = 1, 
i.e., if En(f) = 0(«''). Results of these types are called converse, or inverse theorems. 
S.N. Bernstein pervaded the approximation scene from that time on. The famous 
Bernstein’s inequality, 

IW|sHW." = 0,1,2,.... 

in which r„ denotes a trigonometric polynomial of order n, and which he introduced 
and used in his proof became the model of a key inequality in the investigation of 
converse, or inverse, results on the degree of approximation by best as well as linear 
approximations. Inequalities of this type are called Bernstein type inequalities. 

From the above it follows (see IP. Natanson (1964)) that for it = 1, 2, .... 

Enif) ^ (BiJn’^) coif^'\ m, n = 1, 2, ... ,/G Cjn, 

and that E„(f) = 0(«'**^“^), 0 < a < 1, in turn implies that 5) = 0(5“) and 

moreover, that E„(f) = 0(/i'*^), in turn, implies that 5) = 0(5 In (1/5)). 



18 


The above results are sufficient to conclude that /is infinitely differentiable iff 
liiTin-^ n^E„(f) = 0, for all k. S.N. Bernstein went further to characterize analytic 
functions by their degree of approximation satisfying: limn-+c« (EnCf))’^" < 1, and the 
entire functions by: limn-x» (En(f))^^" = 0. These characterizations were obtained with 
the help of the Fourier series itself. 

The precise characteristic of the structural property of/ that is responsible for 
En(f) = 0(n'^), however, could not get resolved until Anthony Zygmund (1945) 
arrived on the scene and proved that: E„(f) = 0(n^) iff cozif; S) = 0(5), where 

CtbC/i <5) = maxj /,|<5 \\f(x+h)-2f(x)+ f(x-h)\\, 5>0. 

The class of functions satisfying the Zygmund’s condition is known as the Zygmund 
class Z. 

The result of A. Zygmund naturally led to the investigation of the notion of the 
modulus of smoothness of a general order m defined by 

mk(/; 5) = max^|< 5 ||Eo^^ (-1)'"-'' "'Ckfx+kh) ||, (5> 0, 

and the Jackson’s theorem was updated to 

E„(f) 1/n), n = 1, 2, ... ,/e Cm, 

and the Bernstein’s theorem to: 

En(f) = 0(n'“), 0<a<k, 

implies that 

(O,(f;5) = 0i5^). 

The Zygmund’s theorem led to: 

En(f) = 0(n-*) iff (Ok,i(f, 5) = 0(^), k>l, 

the functions satisfying cOk+iif; 5) = 0(<5‘') being said to belong to the generalized 
Zygmund’s class Zk+i- Thus Z is the same as Z 2 . 

The Jackson-Bemstein-Zygmund theorems seemed to round up the chapter on 
the degree of best approximations. However, several mathematicians, mostly related 
with the Russian school, with A.F. Timan as a chief proponent, were not satisfied with 
the orders n'“ in the inverse theorems. The basic question took the shape of when does 
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it happen that E„(/) is precisely of the order Q)\t,(f; 1/n), i.e., under what conditions can 
one assert that there exists a constant 5^ such that 0)^(f; l/n)< Bk E„(f), so that in view 
of the Jackson theorem: En(f) ~ C0k(f, I/n)? In the max-norm case, which is the case we 
are considering in this thesis, a partial answer was provided by S.M. Lozinskii and 
N.K. Bari and S.B. Stechkin (for which see A.F. Timan (1966)) by the statement that 
it happens if/ satisfies 

{f\u)lu^^']du = 0(m,(/;r)), (r 0) . 

This condition, however, is only a sufficient condition, and is not necessary for the 
equivalence E„(f) ~ 1/n) to hold. It is, nevertheless, necessary for the relation: 

Ev"=o (/) = OiE„ if)), (n ^ oo). 

The necessary and sufficient condition for E„(f) ~ cOkif, 1/n), was obtained 
recently by R.K.S. Rathore (1994). His result states that: 

Enif) ~ (Okif, 1/n), in oo) iff (Ok+iif; S) ~ (Okif; S), (5 0), ik > 1). 

For a corresponding development in the degree of approximation by linear 
operators we refer to the texts in approximation theory mentioned earlier, and in 
particular to the work of Z. Ditzian and V. Totic (1987). For the developments in the 
inverse theorems of a general order (pil/n), as compared with the orders in the 
above, both in best approximations and in approximations by linear operators, we 
refer to R.K.S. Rathore (2000) and the bibliography therein. 

An order function (pit) in the thesis is supposed to be associated with a growth 
function ^t) subject to the conditions (see R.K.S. Rathore (1994)) that (p and (ji are 
positive functions on (0, c], such that for all h e (0, 1]: qih) = sup,g(o,c] 
i<pith)(pit))/i(j)it)(pith)) < oo, uih) = sup,e(o,f] (pit)/(pith) < oo, wih) = sup {«(?) : / € [h, 
1] } < oo, and lim/.^ qit) = 0. 

The growth functions 0(r) in the thesis are of the type 0it) = f.ipt = m, or /t = 
2m, m being a positive integer: 1, 2, 3, .... For this 0(O> some examples of appropriate 
(pit) are f, t“log(l/r), r“/(log(l/f))^ 0 < a< pt, any A:-th order modulus of smoothness 
(j)kit) ik < fi) (see G.G. Lorentz (1966)) of any function in L, (1 < ^ < oo), C 2 n etc. 



It may also be remarked that all the above degree of approximation results on 
the best approximation by trigonometric polynomials have been obtained essentially 
by a proper use of appropriate sequences of linear operators themselves that provide a 
so-called “good” trigonometric polynomial approximation (see, e.g., H.S. Shapiro 
(1969)). Thus the history of degree of best approximations is not really separate from 
the history of approximation by linear operators. Such results of our interest, about the 
operators that we have considered in the thesis, are included in the first sections of the 
respective chapters. 


0.4 CONTENTS OF THE THESIS 


The rest of the thesis consists of four chapters: 

Chapter 1 : is a study of the de-blurxing when the blurring operator is modeled 
as the Fej^r sum operator On given by 

(X ” 

(^) = ^ + S & “ « cos kx+b„ sin he ) , 

^ k=i 

where the a* and bk denote the Fourier coefficients of the function / and could be 
expressed in an integral form as: 


<T„(/;x) 


1 


2n7t 


\m 


^ . n{t- x) 
sin 


. t-x 
sin 


dt, n = 0,l,.... 


After a brief history of the operators in Section 1.1, we state some basic results 
required in the proof of our main results about the degree of approximation in the de- 
blurring by 

*=i ) 

in Section 1.2. 

The main degree of approximation results, presented in Section 1.3, of this 
chapter utilize the Peetre’s iiT-functional defined by: 
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where l|/|L = II St"A.(/,^) II. =1 h"/‘"’|L. 

1 

and ||//y denotes the norm of the m-th derivative of the conjugate function of/ 
when m is odd and equals the norm of the ordinary m-th derivative of/if m is even. In 
the latter case, as is well-known, this /f-functional is equivalent to the m-th order 
modulus of smoothness of/. The main results are: 

THEOREM (DIRECT) 

There exists a constant Cm such that 1| 0’„ „f - / |1 ^ C^K^(n~"', /), /e € 2 ^. 
THEOREM (INVERSE) 

I^t the order function (p(t) be associated with the growth function <p{t) = f'. Then: 

pnJ - /II = omyn)) iff Kjrj) = 0((p(t )) , and, 
p.mf-f\\ = o((pil/n) iff Kjr,f) = o((pit))- 

COROLLARY 
Im 0 < Ct <m. Then: 

IK./ -/I = 0(r.-") iffKjrj) = 0 {n,md, 

IK../ - /II = oC"'" ) iff Od" ) . 

Finally, Section 1.4 contains results on numerical simulations with the 
operators <7n,m- We show by an example that a direct inversion does not provide any 
meaningful de-blurring with an actual MRI cross-section. Whereas, on the other hand, 
the iterative de-blurring by the operators cr„,m does provide sharpened images. In the 
next simulation, to obtain the numerical behavior of the error remaining after the de- 
bluiring we take the cross-sections given in the Fig. 0.4 as original cross-sections. 
These we blur by using Fej^r sum blurring and then de-blur by the sequence of de- 
blurring operators Gn,m- The resulting images, the error graphs and the error tables 
have been presented. To assist in understanding the behavior of the error from point to 
point, we also include a graph of the point-wise errors occurring along the horizontal 
midline at the 128* pixel of a cross-section. The graph shows that the errors decrease 
more or less monotonically point wise with the increase in the iteration number. The 
points where the errors are more may be seen to be the points where cross-section is 


less smooth. 
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Chapter 2: studies de-blurring associated with the generalized Jackson 
operators regarded as a blurring model. The generalized Jackson operators are given 
by 


I . n(t-x) 
n sm ^ 


^n.p(f — f f(0 — dl =-^0 + ^Pk''’(^K coskx + bi^ sinAx), 

Ai„ „ ^ t X 2 


for n, p positive integers, with 


''..p = / 


_ sin nr/2 

-U sinr/2 


dt , pr 


where 




= S (-1)^ 

7=0 


2p-l 


, n>l, k>0. 


In Section 2.1 we present some known results on these operators primarily due 
to F. Schurer and F.W. Steutel (1967), R.K.S. Rathore (1974 & 1976), and G.G. 
Lorentz (1966). This is followed by the inclusion of some lemmas necessary to prove 
our main results in Section 2.2. The main results on the de-blurring by the operators 


(-1)' 


,(/;x), n = 1,2,3, ..,(l<n 2 <p 


obtained in Section 2.3, are the following : 

THEOREM (DIRECT) 

For all f € Czn- there holds the error estimate: 

|4"!,/ - /|| s (/;"■' ) . 1 s m s p - 1 , (" a 1). 

where the constant Am.p does not depend on n or f. 

THEOREM (INVERSE) 

Let the order function (p(t) be associated with the growth function p{t) = where 
J :^< p- 1. Then, for f£ C 2 n'. 

\j^:],f-f\ = Omin)) iff C0,Jf;t) = Oi(pil))^and, 
p!!'lf-f\\ = t)((pilfn)) iff (y2„(/;r) = o(<p(r)). 
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COROLLARY 


Let 0 < a < 2m, and Then for f e Czn- 

- /I = 0(n-« ) iff CO,, (/ -,0 = 0(1"), 

iff fflj.(/;r) = o(»“). 

These results are followed by Section 2.4 on numerical simulations wherein 
we consider the cases p = 2, n = 65; p = 3, n = 45; and p = 4, n = 35 of the 
generalized Jackson operators L,,^. These make np-p, the order of the blurring 
polynomial kernel or point spread function (psf) equal 128, 132 and 136, respectively. 
A graph of pint-wise errors along the horizontal midline at 128^ pixel corresponding 
to the same cross-section as considered in the previous chapter confirms the more or 
less monotone decrease of the point-wise errors with increase in iteration number, 
also in the case of the operators of this chapter. 

Chapter 3: considers De La Vallee-Poussin integrals 


V„ = V„if-,x) = 


(2n)!! 


2;r(2«-l)!!4 




2n^-X . 

cos dt, 

2 


as the model of the de-blurring and is a study of the rate of convergence of the de- 


blurring sequence defined as: 




After presenting some well-known results about these operators in Section 3.1, 
in Section 3.2 we obtain some basic lemmas of interest in the proofs of our main 
results. Section 3.3 contains the main results consisting of a direct estimate and an 
inverse result on the corresponding de-blurring operators. These are: 


THEOREM (DIRECT) 

For all f € C 2 m there holds the error estimate: 

- /|| < (/; ) , 1 < m , (« > 1 ). 

w'here the constant Am does not depend on n orf. 


THEOREM (INVERSE) 

Let (pit) be an order function associated with the growth junction (pit) = Then, for 
f € Cin' 

|K‘'"V - f\\ = Oi(pi\l4n)) iff 
jjvW/-/|i=:o((p(l/Vn)) iff <02m(/i5) = o(<)!>(<5)). 
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COROLLARY 

Ler 0 < a < 2m andf€: C 2 n. Then: - /|| = 0(n'^^ ^ ) iff co^Jf ;5) =0(8“ ), 

and, ||K""V-/|| = o(/ 2 -“'^) iff co,^(f:8) = o(d^). 

The numerical simulations are presented in Section 3.4, where in view of the 
fact that the Vallee-Poussin integrals induce a considerable blurring, we utilize their 
filter coefficients coiresponding to the case n = 16x128 whose square root is 
approximately 45, which is the n for the Jackson operators considered with the case p 
= 3 in the previous chapter. Numerical simulations show the behavior of the de- 
blurring errors similar to the case of Jackson operator for p = 3 and n = 45. As in the 
previous chapters, a graph of point-wise signed errors, corresponding to the same 
cross-section as considered in the previous chapters, confirms the more or less 
monotone decrease of the point wise errors with increase in iteration number in the 
operators } . 

Chapter 4: entitled “Applications to MRI diagnosis” consists of numerical 
simulations done in order to obtain an insight in to the applicability of the de-blurring 
in the actual practice where the knowledge of the system blurring operator could be 
obtained only indirectly. In the first section we analyze the filter coefficients of the 
operators considered earlier with the purpose of grading the blurrings induced by 
them. The second section presents numerical studies on cross-filtering in which each 
of the three operators of previous chapters is used in blurring the images. Each blurred 
images is subsequently iteratively de-blurred by all of the associated de-blurring 
operators and the resulting images are cross analyzed. This methodology, with 
appropriate models of blurrings, could be used as an empirical basis for estimating 
system blurring on a particular MRI machine. The experiment also indicates that 
neighboring models perform more or less equally well. 

In the third section we examine the behavior of errors in de-blurring on images 
reconstructed with varying data sizes to see the effect of de-blurring on machine scans 
at different resolutions. Starting with a 512x512 machine data we generate the 
associated images of sizes 64x64, 128x128, 256x256 and 512x512. These are images 
that would be obtained on the machine if the reconstructions of the respective sizes 
were to be acquired with the data so gathered. These are next bluired by the Fejer 



filter associated with an image of size 512x512 by operations in the frequency 
domain.. This is followed by the iterative de-blurring of the same. The errors in 
reconstructing the original image from these (on the 512x512 grid using zero padding 
in the frequency domain) have been calculated. The results show a decrease in the 
errors with the increase of the data size. The de-blurring errors are also calculated by 
regarding the corresponding 64x64, 128x128, 256x256 images themselves as original. 
Tlie results of direct filter inversion (by the same 512x512 Fejdr filter) and iterative 
de-bluiTings performed on the 512x512 machine image are also presented. 

The fourth section in this chapter examines the use of iterative de-blurring to 
obtain Ti, Tt and proton density weighted images from which the calculation of the 
conesponding Tj, T 2 and density values is expected to give more accurate results. We 
start with synthetic Tj, T 2 and density maps used as inputs to the synthesis of Tj, T 2 
and density weighted spin echo MRI images. These images are blurred using the Fejdr 
sums. Subsequently the tissue parameters have been calculated from the blunred 
images and the corresponding de-blurred images. The tabulated results display the 
decrease of the errors in the tissue parameter calculations with the successive de- 
blurrings. 

The chapter concludes by presenting the graphs of the discrete point spread 
functions (the result of applying an operator on an array of data having the central 
entry 1 and the rest zero) of the blurring and the de-blurring operators considered in 
the thesis and by relating their shapes to their numerically observed approximation 
and de-blurring properties. 



CHAPTER 1 : DE-BLURRING BY FEJER SUMS 


In this chapter we study the degree of approximation in the de-blurring by the 
Fejdr sum filters 


m 




m 



In the first section entitled “Fejer Sums” we collect some material about the 
operators {cr„ } , the second section contains some basic results needed in the 
subsequent analysis, and the third section contains our main results on the rate of 
convergence of the sequence The final section on the simulation results 

includes numerical illustrations of the de-blurring achieved on MRI cross-sections. 


1.1. FEJER SUMS 


The complex Fourier coefficients of / e C* are defined as: 

2;r4 

and the real Fourier coefficients are given by: 

=a^if) = — jfit)coskt dt, b,, = bf^(f) = — jf(t)sinktdi, k = 0, 1 , 2 ,..., 

^ -n ^ -n 

with the Fourier series of / as: 

/- 2 =-^^^ + XK(/)cosfcx + Z7i(/)sinfcr] = ]^Ai(/,x). 

(tus ^ 0 

The partial sums of the Fourier series of/are given by 
4 (x) = S„(f;x) = 5 / (x) = ^f(k)e‘^ = ^ (a, (/) cos kx + b,(f) sin kx ) . 

*=-n ^ i=l 

They are expressed as a convolution of / with the Dirichlet kernel as follows: 
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Sin 


(2n + l)(r-jc) 


. t-x 
sin 


-dt, « = 0 , 1 ,..., 


The Fejer sums, or the (C, l)-sums for a function/ e C2n are defined by: 
a,U) = a' (x)= W■^S,(x} + ... + S,,,(x) , £9. + y (i _ 

n '? ^ 

and could be expressed as: 




1 n, oAAi 

= - — f/(r) ^ 

Inn-’ . t-x 

sin 


. n(t - jc) 
sin- 


dt, n = 0,l,..., 


The celebrated result of L. Fejer is as follows: 

1.1.1. THEOREM (L FEJER (1904)) 

If f(x)G C2„ , the sums <T„ (x) converge uniformly along the entire axis toward 
fix): 

lima„(x) = /(x). 

n— >00 


The following well known corollary of the above convergence is due to the 
fact that (C, 1) is a regular summability method: 

1.1.2. COROLLARY 

If the Fourier series of a function fix)E converges at a point Xg then its sum is 

equal to fix^). 

The class Lip m «, 0 < a< 1, defined as: 

LipA/a= { /e C2n: l/(x)-y(y)|<Mlx-y|“ x,yGlR}, 
where IR denotes the real line. 
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The order of approximation by Fejdr sums for the class Lip a,0 < oc < 1 is 
given by: 

1.1.3. THEOREM (S.N. BERNSTEIN (1912)) 

Let f G . Then, for f G Lipj^^C^ , 0< cr < 1 .* |(j„ (jc) - /(x)| <C^M ln“ , x: G iR, the 
factor depending only on a ; and, for f G Lip^ 1 : 
la„(x:) - /(x)| < AM (ln«)/n , x e ^ , n > 1, A being an absolute constant. 


The “In n” term in the above result is a must, as follows from: 


1.1.4. THEOREM (I.P. NATANSON (1964)) 

Suppose that the function f(x)e ^2n has at the point Xq ^ finite derivative 
f'i-Xfy) on the right and a finite derivative flix^) on the left. Then 

r f « ^ Nil flM-f'M 

hm j - — K (^ 0 ) - /(-^o ) J = • 

ln« TT 


In fact, much more precisely, there holds: 

1.1.5. THEOREM (S.M. NIKOLSKII (1940)) 

If A„ia)= sup {max|cr„ (x) - /(x)|} , then: 

feUp,a ^ 

lim„_„ (n / In n)A„ (1) = 2 /;?r , 

and, 

lim„_„ A„ (a) = 2(;r(l - a))"' r(a) sin(a;r / 2) , (0 < a < 1). 

We may mention two other direct error estimates (R.A. DeVore (1972), p.35) 
about the Fejdr sums ct„: 

II f -a fif) |1< (l + ;r)m(/,n-''^), f{x)e Q, 


and 


II f-aXf) P II /'II n-'H, /(Ar)s C»>. 
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where Cj"' denotes the space of 27r-periodic, w -times continuously differentiable 
functions. 


1.2. SOME BASIC RESULTS 


The conjugate series of / is defined by: 


7 ~ ^(-Z?i(/)cosfa:-t-£z^(/)sinfcx) 

1 


If /eC,. the conjugate series of / is the Fourier series of a function / e L^[-n, 7i '\ , 


called the conjugate junction of / . Let H denote the mapping f f , i.e., H f- f . 

In the analysis of de-blurring by the Fejer sums we shall use the saturation 
results about the typical means of the Fourier series. These typical means of order a > 
0, which for 0 < a <2 constitute positive operators, are defined by (R.A. DeVore 
(1972), p. 66): 


Z 1 


1 - 








The m-th iterate cr„ ^ is related to the typical means of the Fourier series by 
a„ m - is easily verifiable by the actions of the two sides on the exponential 

, k =0,±1,.... The following saturation result about{/?“} is due to S. Aljancic 

(1961) (see also G. Sunouchi and C. Watari (1958-59), cf., R.A. DeVore (1972), Th. 
3.5, p. 66). In the sequel we use it to obtain our direct and inverse results about ^ . 


1 .2.1 . THEOREM (S. ALJANCIC (1 961 )) 

The sequence {r“ } is saturated with order n"“ and saturation class 

i.e., \\R“f-f\\=oin~°) iff f is a constant, and I1 /?“/-/ 1|= G(n"“) iff 
feS{R:). 
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The saturation of the sequence of Fej^r sums, due to G. Alexits (1941), 
corresponds to the case a = 1 in above. In view of the saturation of {i?" }, the space of 

zero mean functions / g SiR^) is complete with respect to the norm 

ML =i/iL 

1 

In the sequel, unless otherwise indicated, the norm under consideration is the 
sup norm defined by 

||/||=supl/(x)|, / G Cz„. 

jr 

Certain classes of functions related with our analysis are as follows: 
M denotes the set of functions for which 

|/‘'^||< A/,,0<z<r, 
and 

and 

The class denotes the set of all functions / with |/^''/x)| < M a.e. 

and denotes the union of all (M ) . 

The classes and W‘'’""(l)are same. 

1.2.2. THEOREM (R.A. DEVORE (1972), TH. 1.9, P.16) 

The class of f e Czn ^'hich satisjy^k'A^(f,x)~geL„and ||gL<Afw 
W‘'‘‘/1,M) when r is even and = :(/^"'” ) e Lzp(l,M)} when r is 

odd. 
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The following result is used in the sequel to prove the inverse theorems of a 
general order cp corresponding to the degree of approximation of de-blurring by all the 
filters considered in the thesis: 

1.2.3. Oh-oh-LEMMA (R.K.S. RATHORE (1994)) 

Let c be a positive number, K(l) a non-decreasing positive function and 5(t) a positive 
function on (0, c]. Let (p and <pbe positive functions on (0, c], such that for all h G (0, 
1]: qili) = sup,g( 0 ,,] {4)(Lh)(p{t))l{^t)(p{th)) < oo, u{h) = sup;g(o,c] (p{t)l(p(th) < w{h) = 
sup lu(t) : t E [/i, 1]} < and hm,_,o qif) = 0. Let t„ G (0,c], « > 1, be such that 0 < 

a < tn+i/tn < f3 < I, where a and f are independent of n. If for some constant M, 
K{tj) < M^{t,,) + (p{tj)K(t„)/<p(t^)\, for n > 1 and n < j, then as t 0, 

(\)Kit) = Oi(p(t)) if8{t) = Omt)),(n)K{t) = o{(pit)) if 6(0 = o((p(t)). 

The following two results are also to be used in our subsequent analysis: 

1.2.4. LEMMA 

IfT is a trigonometric polynomial of order n, || T„ ||< 2n |1 T„ ||. 

Proof : Since T = crj + f'ln, the result follows from 1| cr^ |1 = 1 . 

1.2.5. LEMMA 


Proof: From (/)'- sinfcc-i-a* coskx) , (/') = (/)', the result for w = 1 . 

Using induction, assuming the result for m, thus 

^ ^ ((/^'"^ ))' = (( 7 )^'"^)' = ( 7 )''""’\ 

completing the proof. 

We will also need the well-known S.N. Bernstein’s inequality (see e.g., I.P. 
Natanson (1964), p. 90) mentioned in Chapter 0: 



1.2.6. LEMMA (BERNSTEIN’S INEQUALITY) 

IfT is a trigonometric polynomial of order n, || || < n j| || . 


In this chapter, we shall take <p{t) = r", m a positive integer. For this choice of 
(/</), some examples of appropriate qKf) are f, /“log(l/r), r“/(log(l/r))^ 0< a<m, any 
A'-th order modulus of smoothness <D;(r) (A < m) of any function in L, (1 < g < ©o), C 2 n 
etc. In the sequel, the statement “<^0 is associated with means that the order 
function (p(r) and the growth function (pit) are positive functions on (0, c], for some c 
such that for all h e (0, 1]: q{h) = sup,6(o,c] i<lf th) (pit))/ ((p(t)(p(th)) < «>, u(h) = sup,g(o,c] 
(p('t)/(p('th) < oo, H’(/z) = sup {«(/) : r e [A, 1 ] } < CO, and Iimr_>o git) = 0. 


1.3. DEGREE OF DE-BLURRING BY FEJER SUMS 


The Peetre’s K-functionals associated with the typical means of Fourier series 
of integral order are defined as: 

“ .(1/ - *1 + '1«1. ) « J"f (1/ - s|| + - Iff '-’I ) . -« > 1. 


seS(«:) 


g^S(K) 


where H :f —^f denotes the conjugation operation defined before. 


A Jackson-type direct error estimate on the degree of approximation by {cr„ } 
can be given using the K-functional as follows: 

.3.1. THEOREM (DIRECT) 

There exists a constant Cm such that || cr„ ^/ - / || ^ /)» / ^ ^ 2 n- 


Proof : In view of the saturation Theorem 1.2.1, there is a constant Bm such that for 
any ge S(R^), 


(n-r K-i8-S 


<B, 



Hence, „ being bounded, for a constant Cm- 



The result now follows by taking infimum with respect to g £ S(R " ) . 


1.3.2. COROLLARY 

Letf^ C 2 n-If K„(t"',f) = 0((p(t)), then ||a„„/-/|| = C>(<)!)(l/n)), 
and, 

= then ||cr„ „/-/|| = 0 (^( 1 / n)) . 

1.3.3. THEOREM (INVERSE) 

(pit) be an order function associated with the growth function (pit) = f. Then, for 
/e C2„: 

IK.. / -/II = Oi(pil/n)) iff Kjr,f) = Oi(pit)) 

and 

K./-/i|=.«i!a/.) iff Kjr,f)=o(m). 

Proof : By the Lemmas 1.2.4, 1.2.5, and 1.2.6, there exist constants A and B, 
independent of/ and g such that 

Taking supremum over g, 

K. «"./)< la,„f - /|i + 0((r .)" X. 

Hence, for a constant M independent of f, 

Kjrj) < M [\o„J-f\Htri)"KpJ""<f)l 
By the Oh-oh Lemma 1.2.3, K-m/ “/|l ” 0{(pil/n))=^ ,f) = 0{(pit)), and 

K.«-^ “ A\ = o((pil/n))=^ K^(t'”j) = o{(pit)). 

Since the converse is already covered by the Corollary 1.3.2 before, the proof is 
complete. 


As a particular case, taking (pit) = r”, we have: 
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1.3.4. COROLLARY 

Letfe C 2 k, and 0 < a<m. Then: 

Pn . / - /II = 0 (n-“ )iffKjr,f) = Oif ) , 

and, 

- /i =<>(«■“)# a:. «■,/) = 0(1“ ) . 

Note that, as a special case of Theorem 1.2.1, with c denoting a constant, the 
saturation of {a„,m} could be summarized as: 

= L„ , & lltT„^/-/ 1| = o(n"")zj9"/=c. 

Also, in view of the equivalence of ATam and o> 2 m, mentioned in Chapter 0, 
some of the above results for even order iterates may be rephrased as corollaries 
below: 

1.3.5. COROLLARY 

There exists a constant Bm such that |{ (7„ 2 m / “ / II - ^m^im (/J ” ' ) > / ^ 

1.3.6. COROLLARY 

Let (p(t) be an order function associated with the growth function ^t) = Z™. Then, for 
/e C2n- 

||^n 2 m/ - /|| = 0{(pilln)) iff ( 02 „{f-,t) = O((J3(0) 

and 

||<^n.2m / - /|| = I ")) W (f'd) = 0((p(t)) . 

1.3.7. COROLLARY 

Letf e C 2 n, and 0< a<2m. Then: 

\K2J - /II = 0(n -^ ) iff (/;0 = ) > 

and, 

Pn 2«/ - /|1 = "2m (/^ 0 = 0(t“ ) . 



1.4. SIMULATION RESULTS 
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This section contains results on numerical experiments with the de-blurring 
induced by the Fejer operators . The operations of blurring, de-blurring and direct 
inversion are done in both the phase and the frequency encoding directions utilizing 
the two dimensional Cartesian methodology described in Chapter 0. Here: (i) we 
compare the iterative de-blurring with a direct inversion, (ii) compute the 
reconstructions, error tables, and the error graphs for a supposedly known set of 
weighted MRI cross-sections, and (iii) obtain graphs of the point-wise signed errors 
along a straight line through the middle of the cross-section, superimposed on a crop 
of the cross-section. 

1.4.1 ITERATIVE DE-BLURRING COMPARED WITH DIRECT INVERSION 

In the first experiment we verify the effectiveness of the de-bluning scheme 
^n.m compared with the de-blurring based on a direct inversion of the Fejdr filter. 
To start with a model that closely mimics the real life situation, an actual machine 
image output called “Original” is taken to be the example of a blurred image. The 
assumed blurring model is taken to be the Fejdr filter corresponding to n = 128. The 
images associated with the original, four of the iterates cr„^, named “Iteration 1”, 
“Iteration 2”, “Iteration 3”, and “Iteration 4” that correspond to m = 2, 3, 4, and 5 in 
a and the result of the inverse filtering, called “Direct filter inversion” have been 
displayed in Fig. 1.1. 

It is clear from a visual comparison of the images in Fig. 1.1 that the inverse 
Fej^r filter produces severe artifacts, whereas the iterative de-blurrings with „ 
produce much sharper images. The clarity and effectiveness after using the cr„ „, in 
the separation of the gray and the white matter, as well as in delineating the fine 
structure of the entire topography of the regions containing CSF, is amply clear from 
the figure. 



Original 


Iteration 1 



Iteration 4 Direct filter inversion 


Fig. 1.1 ; De-Blurring by Fej6r Suras and Direct Filter Inversion 








1 .4.2 NUMERICAL ERRORS IN THE ITERATIVE DE-BLURRING BY a 
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n m 


The second experiment studies a quantitative behavior of the errors in the 
iterative de-blurrings with „ . To start with a model that closely mimics the real life 
situations, machine outputs associated with images of different tissue contrasts (each 
of which is called an “Original”) are taken to be the examples of actual reference 
cross-sections. Thus, the cross-sections under consideration are assumed to be the 
same as the Ti, T 2 and proton density-weighted two dimensional SE images given in 
Fig. 0.4. The associated gray-levels constitute the functions / for the “Original” 
cross-sections in Fig. 1.2, Fig. 1.3, and Fig. 1.4, respectively. 

These functions / representing the cross sections are blurred by the Fejer 
filter (T„ , with n = 128, in both the frequency and the phase encoding directions, to 
obtain the “Blurred” images shown as the second images in each of the figures 1.2-4. 
The images that result on applying the de-blurring operators „ on these blurred 
images are named “Iteration 1”, “Iteration 4”, “Iteration 7”, and “Iteration 10” and are 
presented in the second and the third image rows in the figures. These reconstmcted 
images correspond to the values m = 2, 5, 8, and 11, respectively in (T^„f . 

A visual comparison of the “Original” images with the “Blurred” images 
shows the nature of blur in the latter images due to the Fejer operator ct„ . The 
progressive return to the sharpness of the original cross-sections through the use of 
^nm seen from the images corresponding to the iterates “Iteration 1” through 

“Iteration 10”, given m the Figs: 1.2- 1.4. 

The Tables 1. 1-1.3 contain the corresponding numerical values of the 
associated errors. The same are displayed in the Graphs 1. 1-1.3 . The percentage 
errors have been, calculated by using the formula: %error = 100x|loriginal- 
reconstructionI|/|loriginall|. where the norms under consideration are the Li, L 2 , and 
L«, norms computed over the entire set of 256x256 pixels in the images. 
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Original 


Blurred 



Iteration 1 


Iteration 7 


Iteration 4 



Iteration 10 



Fig* 1.2 : De-Blurring by Fej6r Sums (IpWeighted Image) 
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Iteration 1 Iteration 4 



Iteration 7 Iteration 10 


Fig. 1.3 : De-Blurring by Fej6r Sums (Ta-Weighted Image) 
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Original 



Iteration 1 



Iteration 7 



Blurred 



Iteration 4 


Iteration 10 



Fig. 1.4 : De-Blurring by Fej6r Sums (Density-Weighted Image) 








Li Error 

L 2 Error 

ho. Error 

blurred 

11 261830 

14 123094 

36 462578 

Ai 

7.380088 

8.570640 

21 222672 

A 2 

6 174041 

6.765291 

16 828613 

A, 

5 501431 

5 805688 

14 031715 

A 4 

5 037621 

5 180599 

12.076420 

As 

4 689713 

4 731967 

10 619979 

A6 

4.416478 

4 390537 

9.486083 

A 7 

4.194215 

4 119931 

8.574721 

As 

4.008613 

3 898783 

7.824581 

Ag 

3.849524 

3 713710 

7.195639 



3 711652 

3.555907 

6.736442 


Table LI : %-Errors in De-BIumng by Fej^r Sums (TrWeighted Image) 



Li Error 

L 2 Error 

Lo. Error 

blurred 

9 862069 

13.041315 

48.783108 

A, 

7.433221 

9.534259 

31.627174 

A2 

6.494610 

8.031424 

24.256310 



5 909182 

7.100663 

19.793056 

A4 

5.488180 

6.446733 

16.759327 

_A 5 

5.162197 

5.959298 

14.542024 

A 6 

4 899826 

5.580925 

12.840414 

A7 

4.682350 

5.274752 

11.488226 

Ag 

4 498209 

5.019303 

10.385134 

A9 

4.338958 

4.801332 

9.466558 

Aio 

4.199403 

4.612151 

8.688783 


Table L2 : %-Errors in De-BIurring by Fejer Sums (T 2 - Weighted Image) 



Li Error 

L 2 Error 

L „ Error 

blurred 

7.287275 

9.150213 

32.423275 

Ai 

4.688092 

5.572545 

23 089937 

A2 

3.913486 

4 438324 

18.803730 

A, 

3.477801 

3 828713 

15.780841 

A4 

3.180496 

3.426594 

13.955763 

As 

2 958243 

3.135305 

12.585967 

A6 

2.783412 

2.912204 

11.490121 

A7 

2.641571 

2.734616 

10.578156 

Ag 

2.523775 

2.589091 

9.799665 

A, 

2.423231 

2.467083 

9.123506 

Ajl! 

2.335815 

2.362873 

8.528898 


Table 13 : %-Errors in De-Blurring by Fejdr Sums (Density-Weighted Image) 
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Graph 1.1 : %-Error in De-Blurring by Fej6r Sums ( Ti-Weighted Image) 
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Graph 1.2 : %-EiTor in De-Blurring by Fej6r Sums (Ta-Weighted Image) 
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Graph 13 : %-EnoT in De-BIuning by Fej6r Sums ( Density-Weighted Image) 



A look at the Table 1.1 reveals that the %-Li-eiTors suddenly drop from 1 1 26 in 
the case of blurred image to its about 66% value 7.38 right in the case of the first iterate. 
The subsequent decrease of the errors is rather slow but monotone. The final error 3.71 in 
the 10th iterate is about 33% of the initial error. The corresponding percentage drops in 
L 2 -errors starting at the value 14.12, are from 60% to 25%, and for the Lo= -error, starting 
at 36.46, the decrease is from 58% to 18%. 

In the Table 1 .2, giving the errors in the case of T 2 -weighted image, the respective 
figures in the Li-case are: 9.86, 75%, and 43%; in the L 2 -case: 13.04, 73%, and 35%; and 
in the L„-case: 48.78, 64%, and 18%. 

Similarly, from the Table 1.3, giving the errors in the case of proton density 
weighted image, these are - in the Lpcase: 7.28, 64%, and 32%; in the L 2 -case: 9.15, 
61%, and 26%; and in the L.«,-case: 32.42, 71%, and 26%. 

The above values are summarized in the following Table 1.4. On the basis of the 
initial errors it gets reiterated that the density weighted images are the smoothest, 
followed by the T 2 weighted and the T] weighted images, which are approximately of the 
same smoothness. The large entry 48.78 in the T 2 -case is probably due to the exceptional 
behavior of the L»-norm in discrete situations (as a single pixel may alter the norm 
drastically and also the numerator and denominator functions may possibly attain their 
maximum values at different pixels). With respect to the Li and L 2 norms (which, 
essentially, are averages), the perfonnance of the de-bluiring, on the whole, is similar for 
the density and the Ti weighted images, followed by the T 2 weighted images. For the Lo, 
norm, however, the gradation of images according to the performance of the de-blurring 
is Ti, T 2 , and density. These gradations, however, are not very pronounced. 
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Table 1.4 : Effectiveness of De-Blurring of Ti, T2, & Density Weighted Images by Fejdr Sums 
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1 .4.3 SIGNED ERRORS IN THE ITERATIVE DE-BLURRINGS BY a„ „ 

n,m 

To get some information on the point-wise behavior of the convergence of the 
iterative de-blurrings by cr„^/, in the following strip graphs(Graph 1.4.(a)-(b)) the 
point- wise signed errors f-CF^.mf in the reconstructions cr^ „/ of / are shown in the 
form of strip-graphs for the first T i-weighted cross-section. These errors were calculated 
for the horizontal 128* mid-row and for all the 256- pixels on it. The unit for the graphs 
is single gray level when the original image / is a 256-gray level image. 

The signed errors for the blurred image are shown in the first strip-graph. These 
errors for the de-blurred cases of the first ten iterates occur in the next ten strip-graphs. 
The last strip-graph is a superposition of all the signed errors of the ten iterates as well as 
the blurred image. The signed differences have been restricted to the range -15 to 25. 
Values crossing this range are not shown. For a better view of the image has been 
compressed in ^-direction. 

These strip-graphs show a more or less monotone behavior of the point-wise 
errors with increasing m. The point-wise errors f -a „^mf ^ along the mid line which is 
shown blackened on the background crop of the image, may be observed to be relatively 
small for the pixels where the cross-section is relatively more smooth and large where the 
cross section is less smooth. The patterns of point-wise signed errors for different m’s are 
more or less similar and the magnitudes of the errors keep on reducing as m increases. 
One may also notice a sudden surge in the magnitude of the point- wise errors if the pixels 
under consideration fall near the boundary of an anatomical structure. In the case of the 
first and the last strip-graphs, some values have actually gone outside the frame and so 
are not visible (due essentially to a limitation of the graphic routine in depicting and 
joining points associated with abrupt and large changes). Such values are associated with 
the sharp peaks occurring at the right. 
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Graph 1.4. (a) : Signed De-Blurring Errors ( / - (T„ „ /, w = 1, 2, 3, 4, 5, 6) 
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Graph 1.4. (b) : Signed De-Blumng Errors ( f - cr„ ,„/, m = 7, 8, 9, 10 ,1 1) & superposition of all 1 1 

above. 




CHAPTER 2: DE-BLURRING BY GENERALIZED 

JACKSON OPERATORS 


This chapter is a study of the rate of convergence of the de-blurring achieved 
by the operators: 

« = 1,2,3,.-, (1 < m < p -1). 

where the generalized Jackson operators {L^ p] are defined as: 


1 ^ 
1 r 


sm- 


\p-n sin — 


dt = -^u:o cos fcc -1- sin fcc) , 


for n, p positive integers and jrs R, the real line, where 




'J sin;/2 I Hi- 


where 


Mr'"=?(-iy 

;=0 


(2p'^np + p-nj-k-\\ 


\ ^ A 


2p-l 


, n > 1, ^ > 0 , 


are the coefficients in the expansion 


f sin nil2\^'’ 


sin r / 2 


np-p 


+2 2,^1''’”^ cask. 


4=1 


The first section entitled “Generalized Jackson Operators” is an introduction 
about the operators the second section contains some basic results needed in 

the subsequent analysis, and the third section contains our main results on the rate of 
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convergence of the sequence ',x)} . The final section on the simulation results 

includes numerical illustrations of the de-blurring achieved on MRI cross-sections 
using {L„ p}. 

2.1. GENERALIZED JACKSON OPERATORS 


The Jackson operators U„ defined by 

3 






Jit) 


. n(t - x) 


-r'l V 


Sin 


. t- X 

sin 


dt. 


were introduced by D. Jackson to prove his celebrated classical Jackson theorem: 


E„if)<C(0(fMn)jG C 2 n, 

where E„{f) denotes the error in the best approximation of / by trigonometric 
polynomials of order n, and co(f;S) denotes the modulus of continuity of /of order 1. 
The generalizations /? > 2, as defined above, have subsequently been considered 
by many researchers, including G.P. Safranova (1950), I.M. Petrov (1958, 1960), P.P. 
Korovkin (1960), I.M. Petrov (1960), Y. Matsuoka (1960), I.P. Natanson (1964), 
Wang Xing-Hua(l964), F.D. van Niekerk (1965), F. Schurer(1965), G.G. Lorentz 
(1966), F. Schurer and Steutel (1967), R.K.S. Rathore (1974, 1976) and many others. 


G.G. Lorentz (1966) used the modification /„.r of the generalized Jackson 
operators defined by 

*=1 


/ =- 1 - 


d cm i- .. 


sin-^ 

L 2 j 


I ; 


]fix-\-kt)dt 


to prove a higher modulus of smoothness estimate: 

fix)-K.rix)\<M,W,M-J^ n = l,2,..., p>i(r + 3) 


R.K.S. Rathore (1974) considered the linear combinations 
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where = det((l //■'“' and A„,+i(Z.nj,) is the determinant obtained from Am+i 
by replacing the first column elements a,i by L,n^(f ; x), i = 1,2, , m+1, and proved 

that: (x) - f (x) || < (/;n ~‘ ), (k< m + 1), where Ca is an absolute constant 

independent of/and n. 

A proof of the basic convergence of Lnj,(f ; x) is given in F. Schurer and F.W. 
Steutel (1967), which also contains most of the results on Lnjif',x) described below. 


2.1.1. THEOREM 

The sequence [L^ p(f',x)} (n = 1, 2, . . .) converges toJ{x) uniformly on [-71, n]. 

The following modulus of continuity estimate for p = 2 is due to P.P. 
Korovkin and was extended for a general p by F. Schurer and F.W. Steutel (1967): 

2.1.2. THEOREM 

For all f (x) G C 2 „ and all integers p>3, 

+ (- = 1 . 2 , . . .) . 


Interestingly enough, the least constant associated with (Oif; l/n) turns out to 
be 3/2. For p = 2, this result is due to Wang Xing-Hua (1964); for p = 3,4 due to van 
Niekerk (1965) and for p > 5 due to F. Schurer and F.W. Steutel (1967). 


2.1.3. THEOREM 

Let p = 2, 3 Then the smallest constant Cp such that 


L„J-f <C^m(/;t)>« = 1.2,...;/GC 


2n ’ 


is given by Cp = 3/2. 


A Voronovskaja-type asymptotic formula for Ln^(f;x) is given by: 
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2.1.4. THEOREM (F. SCHURER AND F.W. STEUTEL (1967)) 


U f ^ ^ second derivative at the fixed point xG [-7t,K'\, then for p>2 


k,p(f'^x)-f(x) = a 




( 1 A 




+ a 




where: 


a = Up -Dip -\) 




flp 
j J 


ip-j) 


2 P-3 


S(p,p) 


p-\ 


2=0 


^2p^ 

y j J 




Here, the functions 5(v, p) arise in the context of the following: 

2.1.5. LEMMA (G.M. FICHTENHOL2 (1964)) 

If V and p are positive integers, then 

S(v,p) = A = t - r -g " §(-i).p^ \2v - 22)=-- . (v>p). 


2^''(2/t-l)!^ 


J ) 


‘r(sinj:)^'' 

S{v,p + i) = \--^dx = 


(- 1 ) 


V+p+\ v-l 


2-^-f2py.f^ 


Y(-^y 


^2v^ 
j J 


p.v-2ff>^ log(2v-2y), 


{v>p). 


For the class of functions Z = {fe C 2 jc- lA.^+O - ‘2 Ax) +Ax-t)\ ^ 2|r|}, which 
actually is a section of the Zygmund-class, 


b„,p =np sup |L„.^(/)-/|, (p^2), 

fuysz 


may be regarded as a measure for the quality of the approximation of the class Z by 
means of the operators L„ ^ . 

For the case p = 2 I.M. Petrov (1958) has proved that 

lim^2„2 =^^^ = 2.647... 

n-^ ’ TT 


For the quantities bp defined by 

Z2^=limZ2,,, (p = 2,3,...) 

;i-4eo ^ 

the following results are due to F. Sdhurer and F.W. Steutel (1967). 



2.1.6. THEOREM 
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If for the operators L„ ^ the quantities are defined as above, then 

S{p,p-\) 


bp =2p- 


Sip,p) 


2.1.7. COROLLARY 



(p-^oo). 


2.1.8. THEOREM 


min b„ 

(p) '' 


K = 


121og2 


7t 


Similarly, for the class Zo = {/e C 2 n- l/(Jc+0 - 2/(x) +/(x-0| t^) (the notation 

Z 2 here is different from that of Chapter 0) of functions there holds: 


2.1.9. THEOREM (F. SCHURER AND F.W. STEUTEL (1967)) 

With p is defined by 


c„,p=n p sup L„,^(/)-/ 


/(x)s2j 


ip^2), 

^2p^ 


c„ = lime 


n,p 




V J J 


ip-j) 


2p-3 


Sip,p) 


gc-i) 


;=o 


. j J 


(P - j) 


2p-l 


2.1.10. COROLLARY 


Cp~ 3 p, (p-400). 


The following results on the approximation of derivatives by the derivatives of 
Lnj,{f-, x) and a Voronovskaja-type formula in this simultaneous approximation were 
obtained by R.K.S. Rathore (1976): 


5FTTmt 

rrljfrni' 

spo a 
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2.1.11. THEOREM (R.K.S. RATHORE (1976)) 

Let p and r be positive integers such that r < 2/7-2. Let f be a In-periodic integrable 
function. If at a point e R, both the right and the left derivatives f+\x) and 
ff^ (jc) exist, then 

+ fy\x))/2. 

Further, iff is such that the derivative f^'^^(x) exists in an open interval containing an 
inten'al [a, i?] c IR and is continuous at each x e [a, b], then uniformly in x E [a, b]\ 

lim„^.L^Z(f;x) = f^^Hx). 

2.1.12. THEOREM (R.K.S. RATHORE (1976)) 

Let p arid r be positive integers with r < 2/7-4. Let f be a 2n-periodic integrable 
function. If at a point x e IR, the right and the left derivatives (tc) and (x) 
exist, then 

4:^ (/ ; X) - / ‘^\x) = (1 - pI” (x) + f^^^ (X)) 12 + o{\ln\n-^oo. 

Further, iffis such that the derivative f^^'"^^\x) exists in an open interval containing 
an interval [a, i»] C IR and is continuous at each x E [a,b'\, then uniformly inxE {a,b]: 

4:!, {f\x)-f^^\x) = {\- I )f {x) + o{\ln\n-+oo. 

It was also pointed out in R.K.S. Rathore (1976) that by the technique of linear 
combinations of the operators obtained in R.K.S. Rathore (1974), the convergence to 
the value f‘''^^{x) can be substantially accelerated for sufficiently smooth functions. 


2.2. SOME BASIC RESULTS 


In this section we collect some basic results needed to prove our results on the 
degree of approximation in the de-blurring by the operators 
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2.2.1. LEMMA (R.K.S. RATHORE (1974)) 

Let m be a positive number and 2/? > w + 1; then 

where Cm.p is a constant independent of n, and ^ ;o)< C 

2.2.2. LEMMA 

= 0(n'^^), n-^oo, 7 = 1,2, 3, ..., (k<p) . 


Proof . For 7=1, the result follows from the previous Lemma 2.2.1. For 7 > 2, by 
induction, assuming the result up to (/-l)-th iterates of L„^, we have 
t-x . . . 2 mt-s^s-x 


;x) = LCi(^n..(sm 




r;-l/r //• ^-S S-X t-S . S - X .2m ^ \ 

^^n.p(^n,p((sin^-cos-y- + cos-^ 


<Li:;(L„,p((2sin^^+2sin^^)'";5);A-) 


Z^A=0 k 




t-s 


\sy,x) 


<2'"y'" f^'V" 

k p-” 


Ik s-x 


,l((sin );x)L„^((sin 


2 (m-/:) 


t-x 


\x) 




m 


t-x 


< 2'” j;” f ^ b(n'“ )C>(/i"^<'"'^' ) = OCn"^™ ) 

* I ^ j 


as the variables t and s are just the dummy variables, completing the proof. 

2.2.3. LEMMA (G. FREUD (1967)) 


2.2.4. LEMMA 

There exist trigonometric polynomials Tb such that for f e , 

\fit)-fix)-'2t;\P\x)/k\]T,(t-x)\<Alf^^'"^\\sin^^ 
where A is a constant independent of f. 
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Proof . For the functions 0!^'^ given by 

Q%,x) = -A^%,x)/A^°\x,x), 

where A^'^\t, x), ^ = 0, 1, ..., m+1, are the determinants of (m+l)-th order defined by 
A^^\t, j:) = I (fly |, where with [/72] denoting the largest integer not greater than jl2. 



1, 

7=1 

<=• 

sin[y/2]r. 

j is even , 


cos[y72]r. 

j is odd 


1, 

7=1 

sin[772]x. 

j is even. 


cos[y72]x. 

j is odd 


■■ l?i 

/ A: + 1, 


R.K.S. Rathore (1974) showed that x) = {t-x'f' + and moreover that 


where, with hxix) = 0, hx{t) is bounded and is continuous sxt = x. From this by using 
L’Hospital’s rule we find that: \hx{f)\ < CflAjc)] + \f"{x)\ + ... + + |lf ^^^||], for 

some constant C, not depending on/or or. 

Replacing m in above by 2m, and taking Tk(t) = 0), in view of the G. 

Freud’s inequality: |lf Ilf ^ ,2m> k>2, and elementary fact 


that for 0<|/-jr| <;r, sin“(r-A:)/2>;r ‘(/-x)^, the result follows. 


2.2.5. LEMMA 

If T is a trigonometric polynomial, IjL^^pr-r =0{n~^"'),n 


oo ^ 


Proof If fk(x) is one of the functions cos kx, or sin kx, in view of l-p"”'’ = 0{n~^), 
which follows from the asymptotic formula of F. Schurer and F.W. Steutel (1967), 

J^:lfk -fk= -a - pr )'”/*= ), n oo . 

Since any trigonometric polynomial is a combination of these functions and the 
operator maps 1 into 1, the result follows. 
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The modulus (Om of smoothness of order m occurring in the sequel is defined 


by: 


= max 


lA|<5 






f{x + kh) 


the norm being the sup-norm with respect to the variable x. 


The following general result using cOmif, 5), was proved by G. Freud (1959) 
for the case m = 2 and was conjectured by him for arbitrary integer m. In full 
generality it was subsequently proved by G.I. Sunouchi (1964): 


2.2.6. LEMMA (G.I. SUNOUCHI (1964)) 

Let B be a linear operator C 2 n —> C 2 n the restriction of B to the space Cj™’ of 

2n-periodic functions having a continuous 2n-periodic m-th derivative. Then, for an 
arbitrary f e C 2 n cind an arbitrary integer v, we have for all x. 


||5(/)||c = max I B(f; x) |< K„, (||S|| + v ' 


B 


(m) 


)coJf;v-^), 


where the norms of the linear operators B, respectively, B^"'^ are defined as 




B 


(m) 


= sup 






We will also use the following version of Jackson-Bemstein theorem (a 
consequence of Theorem 5.1, in R.K.S. Rathore (2000)): 


2.2.7. THEOREM 

Let E„(f) denote the error in approximating f e Czn by the trigonometric polynomial 
of order n of best approximation. Then: 

E„{f) = 0((p(yn)) iff (0,„{f-,t) = O((p(t)), 

and 

E„if) = o((p(\/n)) iff (0,Jf-,t) = o((p{t)). 

In this theorem, and in the rest of this chapter, (p(t) denotes an order function 
associated with the growth function 0(r) = as in the Oh-oh-Lemma 1.2.3. 
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2.3. DEGREE OF DE-BLURRING BY GENERALIZED 

JACKSON OPERATORS 

Based on the results of the previous section, a direct theorem on the degree of 
approximation of/ by the de-blurring operation iJ/J,/ is as follows: 


2 . 3 . 1 . THEOREM 

For all f e C2n , there holds the error estimate: 

pn.lf - /|1 ^ (/;«■’ ) > 1 < m < p - 1 , (n > 1). 

where the constant Am, p does not depend on n or f. 


Proof . Operating on the estimate 

|/(0-/W-S'r;'[/'"’W/^!]T,(t-x)l<A||/'^™>||sin^'"(r-x)/2, 

obtained in the previous section, we have 


Jt=l 


{ 2 k)\ 


In view of G. Freud’s inequalities / 


(t-i) 


A, + ^(2” - l|/'''">j|o(n-'"’ ). 

, ^ > 1, / e C'2 ^ , therefore 


< 2 it 


||zi;J/-/||<fin-^'"||/^'"’||, \<m<p-\, fecg"'\ 

where B is an absolute constant, independent off and n. 


n,p 


■■ 1, (7 > 1) , we have 


Using G.I. Sunouchi’s result, as 

Talcing v = n , the theorem follows. 
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Note that since the approximations are trigonometric polynomials of 

order np-p, their above order of approximation is the same as that of the trigonometric 
polynomials of best approximation given by the classical Jackson theorem: E„(f) < 
CkOhm(f;yn). 

The estimate of Theorem 2.3.1 and the Theorem 2.2.7 of the previous section 
now implies the following: 


2.3.2. THEOREM (INVERSE) 

/f/G C 27 t , 1 ^ rn < p-1 and (pit) is an order function associated with the growth 
function (l)(t) = t^, then: 

\liflf-f\\ = 0((p(l/n)) iff co,Jf-,t) = Oi(p(t)) 


and, 


liflf ■ 


n.pJ 


f ^oicpilfn)) iff o),Jf-,t) = oi(pit)). 


Proof : The ‘if part’ of the result is a consequence of Theorem 2.3.1 and the ‘only if 
part’ that of Theorem 2.2.7. 

A self contained proof of the ‘only if part’ of the above theorem could be 
based on the Bernstein’s inequality 1| Tj || < « 1| T„ || for a Trigonometric polynomial 
of order n, along the lines of the last theorem (Theorem 3.3.2) in the next chapter. 

Taking <p(t) = in Theorem 2.3.2, we have the following: 


2.3.3. COROLLARY 

Let 0 < a<2m < 2p-2, andf € C 2 n ■ Then: 

k"i/ - /I = 0(n-°) iff co,Jf-,t) = 0(n ■ 


and, 


0-/1 = = Oif) 
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Ongmal Iteration 1 



Iteration 4 Direct Filter Inversion 


Fig. 2.1 : Deblurring by Generalized Jackson operator (p=2) and direct filter inversion 
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Original 


Iteration 1 




Iteration 2 


Iteration 4 


Iteration 3 



Direct Filter Inversion 



Fig. 2.2 : Deblurring by Cjeneralized Jackson operator (p=3) and direct filter inversion 
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Oi iginal 



Iteration 2 




Iteration 3 



Direct filter Inversion 


Fig- 23 : Deblurring by Generalized Jackson operator (p=4) and direct filter inversion 
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2.4.2 NUMERICAL ERRORS IN THE ITERATIVE DE-BLURRINGS BY L'"” 

The second experiment studies the quantitative behavior of the errors in the 
iterative de-blurrings with . The setup is essentially the same as that of Chapter 1 . 
The same originals as before are taken to be the examples of actual reference cross- 
sections, which were the Ti, Tj and proton density- weighted SE images given in Fig. 0.4. 
The associated gray-levels constitute the functions / for the “Original” cross-sections in 
Figs. 2.6-2.8. 

These functions / representing the cross sections are blurred by the Generalized 
Jackson filters L„ ^ (p = 2, 3, 4, with respective n = 65, 45, 35), in both the frequency 
and the phase encoding directions, to obtain the “Blurred” images shown as the second 
images in the Figs. 2. 6-2.8. The images that result on applying the de-blurring operators 
4,”^ on these blurred images are named “Iteration 1”, “Iteration 4”, “Iteration 7”, and 
“Iteration 10” and are presented in the second and the third image rows in the figures. 
These reconstructed images correspond to the values m = 2, 5, 8, and 11, 
respectively. 

A visual comparison of the “Original” images with the “Blurred” images shows 
the nature of blur in the latter images due to the generalized Jackson operators L„ ^ . The 
progressive return to the sharpness of the original cross-sections through the use of 4!^^ 
may be seen ft-om the images corresponding to the iterates “Iteration 1” through 
“Iteration 10”, given in the Figs: 2.6-2.8. 

The Tables 2.5-2.7. contain the corresponding numerical values of the associated 
errors. The same are displayed in the Graphs 2.4-2.6 . The %-errors have been calculated 
by using the formula: error = 100x||original-recoiistruction||/||original||, where the 
norms under consideration are the Li, L 2 , and Leo norms computed over the entire set of 
256x256 pixels in the images. 
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Iteration 1 Iteration 4 



Iteration? Iteration 10 

Fig. 2.4 : De-Blurring by Generalized Jackson operator (p=2) (Ti weighted image ) 
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Iteration 4 
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Fig- 2*5 : De-Blurring by Generalized Jackson operator (p=2) (T 2 weighted image) 
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Iteration I 



Iteration 7 


Iteration 4 



Iteration 10 


Fig* 2*6 : De-Blurring by Generalized Jackson operator (p-2) (Density weighted image) 











Lj Error 

Lj Error 

L Error 

blurred 

11.826023 

15 174283 

39.176308 

Aj 

9.655653 

11 779855 

28.377062 

A2 

8.924964 

10542892 

25 624622 

A3 

8 486915 

9.811369 

23.834412 

A4 

8.168295 

9 295283 

22.535946 

_ As 

7.918212 

8.899959 

21 526028 

^ 

7 714353 

8.582493 

20 702404 

A7 

7.542678 

8.319170 

20 007799 

^ 

7.395671 

8.095625 

19.407467 

A9 

7.266559 

7.902443 

18 878998 

Aio 

7 152111 

7.733020 

18.407181 


Table 2.1 : %-ErTor in De-Blurring by Generalized Jackson operator (p=2) 

(Ti weighted Image) 



Li Error 

Li Error 

L 0 , Error 

blurred 

10.898792 

15.052915 

57.442039 

A, 

9.570345 

13 023810 

45.283543 

42 

9.032103 

12 105190 

39.900150 

A, 

8.706363 

11514212 

36.862049 

. A4 

8.460309 

11.075411 

34.657135 

...As 

8.258093 

10.721021 

32.933746 

A6 

8.087586 

10.421531 

31.522835 

A7 _____ _ 

7.940648 

10.161642 

30.331257 

As 

7.812851 

9.932296 

29.302383 

A, 

7.700998 

9.728980 

28.399206 

Aio 

7.602109 

9.548022 

27.596172 


Table 2.2 . %-Error in De-Blurring by Generalized Jackson operator (p=2) 

(T 2 weighted Image) 



Li Error 

Lz Error 

L „ Error 

blurred 

7.587895 

9.755048 

37.955128 

A, 

6.166692 

i.enm 

32.114063 

Ai 

5.685149 

6.856499 

30.864574 

As 

5.397155 

6.411814 

29.555965 

A4 

5.191302 

6.096819 

28.291136 

As 

5.031132 

5.853706 

27.114025 

As 

4.900396 

5.656725 

26.034811 

A, 

4.789762 

5.492030 

25.050058 

Ag 

4.694689 

5.351151 

24.151577 


4.611548 

5.228565 

23.330147 

Aio 

4.537868 

5.120420 

22.576988 


Table 23 : %-EiTor in Efe-Blurring by Generalized Jackson operator (p=2) 
(Density weighted Image) 
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Iteration 10 


Fig. 2.7 : De-Blurring by Generalized Jackson operator (p=3) (Ti weighted image) 
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Fig. 2.8 : De-Blurring by Generalized Jackson operator (p=3) (T 2 weighted image) 
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Fig. 2.9 : Dc-Blurriiig by Generalized Jackson operator (p=3) (Density weighted image) 









Li Error 

L 2 Error 

L„ Error 

Blurred 
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17.429298 

45.340599 

A, 
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13.792450 

34 178101 

A 2 

10.167757 
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A, 

9.755668 
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A4 
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11.299795 
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As 

9.255418 
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A 6 
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A, 

8.929787 

10.377771 
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24 136005 

A9 
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Ajo 
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9 814838 
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Table 2.4 . %-Error in De-Blumng by Generalized Jackson operator (p=3) 

( Ti weighted Image) 



Li Error 

Li Error 

L„ Error 

blurred 

11.878249 

16.498714 

64.058960 

A, 

10.547678 

14.551584 

53.267841 

Ai 

10.041962 

13.712990 

47.590000 

A, 

9.730365 

13.189909 

43 934147 

A4 

9.513884 

12.812903 

41.970104 

As 

9.354856 

12.520995 

40.485146 

As 

9.229193 

12.284058 1 

39.281612 

P 

A7 

9.123915 

12.084971 

38.271519 

Ag 

9.030957 

11.913160 

37.402035 

* *0 

Ag 

8.946622 

11.759972 

36.639194 

Aia 

8.869592 

11.621891 

35.959953 


Table IS : %-Error in De-Blurring by Generalized Jackson operator (p-3) 

(T 2 weighted Image) 



Li Error 

Li Error 

L„ Error 

Kliirrtf^d 

8.557381 

11.186337 

41.888958 

Ai 

6.994592 

8.855388 

35.729176 



A 7 

6.498222 

8.060233 

33.982243 

- 

A-^ 

6.221765 

7.623755 

33.660534 

. .. 

A 4 

6.032843 

7.328251 

33.222908 

— — 

5.889729 1 

7.105361 

32.732426 

, — 

As 

5.773912 

6.926474 

32.224869 

— 

A 7 

5.677165 

6.777238 

31.719654 

. -- „ ... 

Ag 

5.594154 

6.649395 

31.227001 



Ao 

5.521667 

6 537702 

30 751942 

— 

Ahi- 

5.457621 

6.438704 

30.296625 


Table 2.6 : %-Error in De-Blurring by Generalized Jackson operator (p=3) 
{ Density weighted Image) 
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Fig. 2.10 : De-Blurring by Generalized Jackson operator (p=4) (Ti weighted image) 













Iteration? Iteration 10 

Fig. 2.11 : De-Blurring by Generalized Jackson operator (p=4) (T 2 weighted image) 
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Fig. 2.12 : De-Blurring by Generalized Jackson operator (p=4) (Density weighted image) 








Li Error 

L 2 Error 

Loo Error 

blurred 
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Ax 
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A2 
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13.751358 
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Ai 

10 504880 

13 007795 

30 409939 

A4 

10 226680 

12 514482 

29.230158 

As 

10.018459 

12.148420 

28 390139 

A^ 

9 853121 

11.858305 

27.713596 

A7 

9.716592 

11 618515 
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As 
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11.414445 

26 675184 

Ai 
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1 1 237056 
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Aio 

9.404180 
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Table 2.7 : %-EiTor in De-Bluning by Generalized Jackson operator (p=4) 
(Ti weighted Image) 



Li Error 

L 2 Error 

Loo Error 

blurred 

12.518543 

17.406057 

67.959724 

Ai 

11.127961 

15 453486 

58.020721 

A2 

10.629580 

14.633122 

52.682335 

A3 

10.334711 

14 132834 

49.198338 

A .4 

10.122963 

13.775389 

46.677059 

As 

9.959428 

13.498801 

44.732990 

A4 

9 828265 

13.274008 

43.509285 

At 

9.721177 

13.086461 

42.582291 

As __ 

9.633188 

12.926210 

41 790131 

A, 

9.558999 

12.786831 

41.099354 

Aio 

9.494246 

12.663932 

40 487347 


Table 2.8 : %-EiTor m De-Blurring by Generalized Jackson operator (p=4) 

(T 2 weighted Image) 



Li Error 

Lj Error 

Loo Error 

blurred 

9.244117 

12 198750 

44.153336 

Ai 

7.545388 

9.684868 

38.582577 

At 

7.002904 

8.820059 

35.470993 

A, 

6.717432 

8.360008 

34.856041 

A 4 

6.527386 

8.059769 

34.813709 

As 

6.386542 

7.839794 

34.681816 

A^ 

6.275219 

7.666822 

34.494614 

A 7 

6.183958 

7.524434 

34.274326 

As 

6.106483 

7.403521 

34.035374 

A«) 

6.038725 

7.298472 

33.787163 

Aio 

5.978663 

7.205680 

33.535870 


Table 2.9 : %-EiTor in De-Blurring by Generalized Jackson operator (p=4) 





%-Error %-Error %-Error 



Number of 


Graph 2.7 : %-Error in De-Blurring by 
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Graph 2.8 : %-Error in De-Blurring by Generalized Jackson operator (p=4) (Ta weighted image) 



Graph 2.9 ; %-Error in De-Blurring by Generalized Jackson operator (p=4) (Density weighted image) 
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A look at Table 2.1 reveals that the %-Li-errors drop from 11.82 in the case of 
blurred image to its about 81% value 9.65 in the case of the first iterate. The subsequent 
decrease of the errors is rather slow but monotone. The final error 7.15 in the 10th iterate 
is about 60% of the initial error. The corresponding percentage drops in La-errors starting 
at the value 15.17, are from 77% to 50.9%, and for the Loo -error, starting at 39.17, the 
decrease is from 72% to 46%. 

In Table 2.2, of Ta-weighted image, the respective figures in the Li-case are: 
10.89, 87%, and 69%; in the La-case: 15.05, 86%, and 63%; and in the Loo-case: 57.44, 
78%, and 48%. 

In Table 2.3, of proton density weighted image, these are - in the Li-case: 7.58, 
81%, and 59%; in the La-case: 9.75, 79%, and 58%; and in the Loo-case: 37.95, 84%, and 
59%. These values are summarized in the following Table 2.10. 


1 

Drop in the Lj-eiror from 
the first to the last iterate 

Drop in the Li-error from 
the first to the last iterate 

Drop in the Loo-error from 
the first to the last iterate 

Initial 
% -Error 

First 

Iterate 

Last 

Iterate 

Initial 

%-EiTor 

First 

Iterate 


Initial 

%-Error 

First 

Iterate 

HQHH 

B 

10.89 

87% 

69% 

15.05 

86% 

63% 

57.44 

78% 

48%. 


11.82 

81% 

60% 

15.17 

77% 

50% 

39.17 

72% 

46% 


7.58 

81% 

59% 

9.75 

78% 

52% 

37.95 

84% 

59%. 


Table 2.10: Effectiveness of De-Blurring of Tj, T 2 , and Density Weighted Images (p = 2) 


For p = 3, Table 2.4 reveals that the %-Li-errors drop from 13.30 in the case of 
blurred image to its about 82% value 10.91 in the case of the first iterate. The subsequent 
decrease of the errors continues to be slow but monotone. The final error 8.58 in the 10th 
iterate is about 64% of the initial error. The corresponding percentage drops in L 2 -errors 
starting at the value 17.42, are from 79% to 56%, and for the -error, starting at 45.34, 
the decrease is from 75% to 51%. 

In Table 2.5, of Ta-weighted image, the respective figures in the Li-case are: 
11.87, 88%, and 74%; in the La-case: 16.49, 88%, and 70%; and in the Lc„-case: 64.05, 
83%, and 56%. 
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In Table 2.6, of proton density weighted image, these are - in the Lpcase; 8.55, 
81%, and 63%; in the Li-case: 11.18, 79%, and 57%; and in the Loc-case: 41.88, 85%, and 
72%. These values are summarized in the following Table 2.1 1. 



Drop in the Li-error from 
the first to the last iterate 

Drop in the L 2 -error from 
the first to the last iterate 

Drop in the Loo-error from 
the first to the last iterate 

Initial 

%-Error 

First 

Iterate 

m 

Initial 

%-EiTor 

First 

Iterate 

Last 

Iterate 

Initial 

%-'Error 

First 

Iterate 

B 

m 

13.30 

82% 

64% 

17 42 

79% 

56% 

45.34 

75% 

51% 

a 

11.87 

88% 

74% 

16 49 

88% 

70% 

64.05 

83% 

56%. 


8.55 

81% 

63% 

11.18 

79% 

57% 

41.88 

85% 

72%. 


Table 2.11: Effectiveness of De-Blurring of Ti, T 2 , and Density Weighted Images (p = 3) 


In Table 2.7 the %-Li-errors drop from 14.35 in the case of blurred image to its 
about 81% value 11.75 in the case of the first iterate. The subsequent decrease of the 
errors is slow but monotone. The final error 9.40 in the 10th iterate is about 65% of the 
initial error. The corresponding percentage drops in L 2 -errors starting at the value 19.00, 
are from 79% to 58%, and for the L=o -error, starting at 49.46, the decrease is from 76% to 
52%. 

In Table 2.8, of Tz-weighted image, the respective figures in the Lrcase are: 
12.51, 88%, and 75%; in the La-case: 17.40, 88%, and 72%; and in the L«-case: 67.95, 
85%, and 59%. 

In Table 2..9, of proton density weighted image, these are - in the Li-case: 9.24, 
81%, and 64%; in the La-case: 12.19, 79%, and 59%; and in the L„-case: 44.15, 87%, and 
75%. These values are summarized in the following Table 2.12. 


■ 

Drop in the Li-error from 
the first to the last iterate 

Drop in the L 2 -error from 
the first to the last iterate 

Drop in th 
the first to 

e Loo-erroi 
the last it 

■from 

erate 

1 

Initial 

%-ErTor 


KSI 

Initial 

%-Error 

First 

Iterate 

Last 

Iterate 

Initial 

%-Error 

First 

Iterate 

Last 

Iterate 

g 


81% 

65% 



58% 

49.46 

76% 

52% 

19 


88% 

75% 

17.40 

88% 



72% 

67.95 

85% 

59%. 

m 



64% 



59% 

44.15 

87% 

75%. 


Table 2.12: Effectiveness of De-Blurring of Ti, Tj, and Density Weighted Images (p = 4) 
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2.4.3 SIGNED ERRORS IN THE ITERATIVE DE-BLURRINGS BY 

n,p 


The information on the point-wise behavior of the convergence of the iterative de- 
blurrings by (p = 2, 3, & 4), is given in Graphs. 2.10 (a)-(b)-2. 12 (a)-(b), for p = 2, 
3, & 4, where the point- wise signed errors in the reconstructions of / 

are shown in the form of strip-graphs for the first Ti-weighted cross-section. These errors 
were calculated for the horizontal 128*'' mid-row and for all the 256- pixels on it, as 
before. The unit for the graphs is one gray level when the original is assumed to be a 256- 
gray level image. 

These strip-graphs too, show a more or less monotone behavior of the point-wise 
errors with increasing m. The point-wise errors / - (p = 2, 3 & 4) along the mid 
line, shown blackened, may be observed to be relatively small for the pixels where the 
cross-section is relatively more smooth and large where the cross sections is less smooth. 
The error patterns for different m’s are more or less similar and the magnitudes of errors 
keep on reducing as m increases. One may also notice a sudden surge in the errors as the 
pixels under consideration come near the boundary of an anatomical structure. 

A comparison of these strip graphs shows that the maximum error 
amplitudes occur in the case of p = 4, followed by p = 3 and p = 2. The nature of the 
oscillations is similar in all of these. The greater oscillations occur near the edges of 
anatomical structures, as had been the case with the Fejdr sum operators in Chapter 1. 
However, the oscillations for the generalized Jackson operators seem to diminish much 
more slowly as compared with the Fejdr case. This indeed is because of the relatively 
poor approximation properties of the present operators. For, going by the theoretical error 
estimates, say e.g., for Lip a functions, the respective errors should be roughly 
proportional to a* powers of 1/35, 1/45, 1/65, and 1/128, in the respective cases of the 
Jackson operators with p = 4, 3, 2, and the Fejdr sums. 
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CHAPTER 3 : DE-BLURRING BY THE DE LA 
VALLEE-POUSSIN INTEGRALS 


In this chapter we study the degree of approximation in the de-blurring 
associated with the De la Vallee-Poussin integrals (see I.P. Natanson (1964)): 


K = 


(2n)\\ 


2;r(2n-l)!!4 




cos 


2n 


dt. 


In the first section entitled De la Vallee-Poussin Integrals we collect some 
material about the operators {y„}, the second section contains some basic results 
needed in the subsequent analysis, and the third section contains our main results on 
the rate of convergence of the sequence ;x)} defined as: 


vr'(/;4 = E"., 


The final section on the simulation results includes numerical illustrations of the de- 
blurring achieved on MRI cross-sections. 


3.1. DE LA VALLEE POUSSIN INTEGRALS 


These (see J. De la Valee Poussin (1908)) derive the importance from their 
simplicity, interesting approximation properties and theirs being variation diminishing 
(a result due to I.J. Schoenberg (1959)). Some well known quantitative results about 
them are as follows: 


3.1.1. THEOREM (I. P. NATANSON (1944)) 

If a function fix) e has modulus of continuity co(5) , then the inequality 


^„ix)-fix)\^3co 




holds for all values of x . 
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3.1.2. THEOREM (I.P. NATANSON (1946)) 

If (oc) = supjmaxjV^ (x) - / wherein f (a:) runs through all functions with 

f ^ 

period 2n of class Lip, a, then 


where lim p„ =0. 

n -400 


{/„(a) = -4l=r{^^“^ 




+ 



3.1.3. THEOREM (I. P. NATANSON (1944)) 

If the function f (j:) G € 2 ^ possesses for a value of x a finite second derivative f'ix) , 
then the formula 

K(^) = /(jc)+— +— 

n n 

with lim p„ = 0 holds true for this x . 

jrj— >00 


The following theorem concerns the point-wise convergence of the derivative of De 
la Vall^e-Poussin integrals; 


3.1.4. THEOREM (I. P. NATANSON(1964)) 

If at a point x, f{x) has a finite derivative f'(x), then the relation 

\imv'ix) = f'ix) 

n-^00 

holds true for this value of x . 


3.1.5. THEOREM (P.C. SIKKEMA AND R.K.S. RATHORE (1976)) 
Let 0 < a < 1 , and p be a positive integer. If 

v.am-'Z^--f^(t-xyy,x) 


£(a,p;K,x) = sup- 

/ 


v=0 


s xe R, 


where f runs through all functions with € Lip, a , then with (a) ^ = 0,1, (a + i - 1) , 

± LL2-^<limA£(a,p;K,x)<- . 

(l4-a)-V7t (L + a)pyIn 



3.1.6. COROLLARY 

If p is a positive integer, then 
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lim n 

n — > oo 


(p + l)/2 


E(Up;V„ 


x) 


+ 1 ) 

ip + 1 )! 


The following result generalizes the approximation of first derivative result 
Theorem 3.1.4 of I.P. Natanson (1964) and shows that the De la Valee Poussin 
integrals, in fact, enjoy the pointwise convergence in simultaneous approximation of 
derivative of an arbitrary order: 


3.1 .7. THEOREM (P.C. SIKKEMA AND R.K.S. RATHORE (1976)) 

I^t m be a positive integer. If f s and at a point xG R, f’'’”\x) exists, then 


3.2. SOME BASIC RESULTS 


In this section we derive some results which we use in the proofs of our main 
results in the next section. 


3.2.1. LEMMA 

Letfor k>l, = V„(sin‘* Then: s n r(-y-)n , « ~ . 


Proof . Since, 


n 

\ 
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it follows that 






(n + k) 


n+k+)^ 


n^) 


= 27i-*r(^^), 




r(^^) 


-Jn 


-n 


3.2.2. LEMMA 


If T is a trigonometric polynomial, 


=0(n"”),n^oo. 


Proof : It follows along the lines of the corresponding result about the generalized 
Jackson operators obtained in the previous chapter. 


3.2.3. LEMMA 

Let 0<m<n be integers. Then 

g2m 


dx 


2m 


cos“-i-i = cos^-i^ 5;«W(")tan 

^ j=0 


y t-X 
2 ‘ 


where g^"\n) is a polynomial in n of degree < j . 


Proof : For m = 0 , the result is trivial. Hence assuming it for m and differentiating 
twice 


^2m-f2 


cos 


2n 


t-X 

~T 


;COS 


t-x 


y=0 


[ml /•„) 


(n-;)(n- 7 -|)tan 


’j+2 


t-x 


- {(« - })ij + i) + (n - y + i)y}tan'^ ~ + jU - i) tan'^'^ , 

t X t X 

from which the result is clear, as the coefficients of tan^^^^-y-, tan^^-y-, and 
inside the square brackets are respectively of degree 2, 1, and zero in n . 

2 
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3.2.4. LEMMA (BERNSTEIN TYPE INEQUALITY) 


Ikv,'”'/) 


(2m) 




/€C 2 „, m,n>l. 


where the constant does not depend on n or / . 


Proof : As we have to show that = 6>(n'”)|/|| , in view of the Lemma 3.2.3 we 

just have to verify that 

i — f cos^'""-'’ — sin^-' — dx = 0{n ~^ ) , 0< j<m, 

7r(2n-l)!!i 2 2 

which, however, is a consequence of our Lemma 3.2.1. 


The m-th order modulus of smoothness for fe C 2 n, <5 > 0, /n > 1, is defined 


by 


<y™(/;5)=sup|»i<j 


^ ni^ 


: fix+jh) 


j 


and the m-th order Peetre’s JST-functional for /€ C 2 n, t>0, m > 1, by 


(/;0 = ~ ^1 • 


We may recall the following well-known equivalence (see e.g., Z. Ditzian and 
V. Totik (1987)) between the two: 

Arr,0)m(f; t)< 0, C2n, t > 0, m > 1. 

where A„ and are absolute constants independent of t and/. 

In this chapter, we shall take ^r) = /^, m a positive integer. For this choice of 
some examples of appropriate q){t) are r“ r“log(l/r), r"/(log(l/0)^ 0 < a<2m, 
any ifc-th order modulus of smoothness (Okif) (k < 2m) of any function m (1 < ^ < 
°°), C 2 n etc. As before, the statement '"(pit) is associated with ^r)” means that the 
order function (p(0 and the growth function 0(r) are positive functions on (0, c], for 
some c such that for all h e (0, 1]: g(ft) = sup.e(o,c] {m)mm)(p(th)) < u(h) = 
sup.e(o.ci <P(t)/(p(th) < oo, wQi) = sup {uit ) : r £ [Ii, 1] } < «, and lim,-^) qit) = 0. 
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3.3. DEGREE OF DE-BLURRING BY THE DE LA 
VALLEE-POUSSIN INTEGRALS 

Based on the results of the previous section, a direct theorem on the degree of 
approximation of / by the de-blurring operation is as follows: 

3.3.1. THEOREM (DIRECT) 

For all f E C 2 it , there holds the error estimate: 

where the constant Am does not depend on n or f. 

Proof : Operating on the estimate 

l/(0-/(^)-Sjrt/'*'W/*nr,(r--x)l<A||/''™>||sin"'"(r-x)/2, 
obtained in the Lemma 2.2.4, and using Lemma 3.2. 1 and Lemma 3.2.2, we have 

Jt=l 

From which, in view of the G. Freud’s inequalities given in Lemma 2.2.3, as before 

where B is an absolute constant, independent of / and n. Hence, using the result of 
Lemma 2.2.6. of G.I. Sunouchi (1964), due to ||v/ 1| = 1, {j > 1) , we have 

from which, taking v = , the result follows. 
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For the De la Val€e Poussin integrals, referring to the Oh-oh-Lemma 1.2.3, 
below we shall use the order function ^t) associated with the function 0(t) = 


3.3.2. THEOREM (INVERSE) 

Let (p{t) be an order function associated with the growth function (p^t) 
forfe C 2 n,'- 

\yt^f-f\\ = 0{(p{l/yf^)) iff co,„if;5) = 0(,(pi5)), 


and 



f-f\\=omi/^)) 




i'”. Then, 


Proof: By the Lemma 3.2.4, there exist constants A and B, independent of g 6 , 

and /e C 2 ;r such that 

if,. (/;»’" ) < if,. (/ - ) + if„ (y,i-)(/ + ) 

S|K'"’/ - /|| + e’lAn’lf - «| + B II* “">||] . 

Taking infimum over g e , for a constant M independent of /, 

K,„(f-,t^)<M -41 +((' riTK^^fin")] . 

By the Oh-oh-Lemma 1.2.3, therefore, we have 

K”'/ - 41 = o(«>("‘'"))=> (/;<’” )=o{<p(t)). 

and 


By the equivalence of K 2 m and cozm, and since the direct part of the result is already 
covered in the Theorem 3.3.1, the proof is complete. 


3.3.3. COROLLARY 

Let0<a<2m. Then, forf e € 2 ^: 

lvW/-4| = Ci(n-“'"i iff <o,.(/;5) = 0(5“). 


and 


- /|1 = ) iff (Oj„{f ;5) = o(5“) . 



3.4. SIMULATION RESULTS 
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This section contains numerical simulations with the de-blurring induced by 
the De La Vallee-Poussin integral V„ . As in the previous chapters, the blurring, de- 
blurrings and the direct inversion are done in both phase and frequency encoding 
directions using the two dimensional Cartesian methodology described in Chapter 0. 
As before, we; (i) compare the iterative de-blurring with a direct inversion, (ii) 
compute the reconstructions, error tables, and the error graphs for the supposedly 
known set of weighted MRI cross-sections, and (iii) obtain graphs of the point-wise 
signed errors along the straight line through the middle of the cross-section, 
superimposed on a crop of the cross-section. 

3.4.1 ITERATIVE DE-BLURRING COMPARED WITH DIRECT INVERSION 

The first experiment verifies the effectiveness of the de-blurring scheme Vj"' 
as compared with the de-blurring based on a direct inversion of the De La Vallee- 
Poussin filter. The same actual machine image output, as taken in the previous two 
chapters, called “Original” is taken to be the example of the blurred image. The 
assumed blurring model here is the De La Vall6e-Poussin filter corresponding to n 
= 16x128, out of which only the 128 filter coefficients are utilized in the simulations. 
The images associated with the original, four of the iterates Vj'”’ , named “Iterationl”, 
“Iteration 2”, “Iteration 3”, and “Iteration 4” that correspond to m = 2, 3, 4, and 5 and 
the result of the inverse filtering, called “Direct filter inversion” have been displayed 
in Fig. 3.1. 

It is clear from the visual comparison of the images in Fig. 3.1 that the inverse 
De La Vall6e-Poussin filter produces severe artifacts, whereas the iterative de- 
blurrings with Vj”’ produce much sharper images. Their clarity and effectiveness in 
delineating the fine structure of the entire topography of the regions containing CSF, 
is amply clear from the figure. Note also a better gray and white matter discrimination 
as compared with the case p—3 in iterations 2-4 of chapter 2 (Fig. 2.2). 
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Iteration 4 Direct Filter Inversion 


Fig. 3.1. : De-blurring by De La Vall6e-Poussin and direct filter inversion 
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Original 


Blurred 



Iteration I 


Iteration 4 



Iteration 7 


Fig* 3.2* : De-Blurring by De La ValI6e-Poussin integrals (Ti weighted image) 








Original 


Blurred 



Iierationl Iteration 4 



Fig. 3.3. : De-Blurring by De La ValMe-Poussin integrals ( T 2 weighted image) 
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Iteration 1 Iteration 4 



Fig. 3 . 4 . : De-BIuning by De U Vall^e-Poussin integrals (Density weighted image) 










Li Error 

L 2 Error 

L„ Error 

blurred 

13.518874 

17.757774 

46268749 

jii 

i 1 0(37230 

14.037602 

34 946190 


10.271811 

12 694630 

30083900 

A 3 

9.836020 

11.938289 

28 262102 

A 4 

9.533145 

11.418419 

27 065260 

A 5 

9.M)I29() 

1 1 022747 

26.133268 


9.ri2()67 

10703659 

25.376198 

A 7 

8.951344 

10.436556 

24.741747 

Ah.. ^ . 

SH115(>8 

10.207154 

24.197197 

A^ 

8 61^163 

10 006439 

23 720888 


s .‘tTxo t'i 

9.828348 

23.297813 


Table 3.1 . -brr(tr!i in Dc-BIunring by De La Vall6e-Poussin integral 

(T| weighted image) 



Li Error 

L 2 Error 

L„ Error 

blurred 

12.{K)5098 

16 666893 

64 805298 

Ai , , 

10 635724 

14.673879 

54.058949 

Ai - . , 

10.106552 

13.794727 

48.286686 


9.776073 

13.239385 

44.509380 


9.541142 

12,834176 

42 134560 

Af 

9 t(>4327 

12 516259 

40548008 

Ajfe ! 

224312 

12 255375 

39 251434 

' A 7 

9,107256 

12.034124 

38,155491 

Alt 

{>(M2^)5 

11 841676 

37 206196 

A'* 

8 9T(X,51 

11.670034 

36 368710 

Aju 

K 824517 

11.513990 

35 619312 


Table 3.2 : 'S -Kmirs in Dc-Blurring by De La Vall6e-Poussin integrals 
(T; weighted image) 



Li Error 

Lz Error 

L„ Error 

blurred 

8.6%555 

11.395054 

42.332867 

A 1 

7.093175 

9.006601 

36.220360 

Ai 

6.567901 

8.169998 

33.887886 

.A 

A 4 

6.274099 

7.707120 

33.582581 

6 071949 

7.392573 

33.156796 

Afs 

5.918304 

7.154491 

32.671185 

Ag 

5.793479 

6.962598 

32.161728 

A 

Ag 

5.688051 

6.801755 

31.648882 

5 597083 

^ 6 663292 

31.144051 

.AJla 

Aj, 

5.517218 

6.541728 

30.653318 

A|» 

5.446238 

6.433475 

30.179653 


Table 3J. : %-Errors in De-BIurring by De U Vall^e-Poussin integrals 
(Density weighted image) 






• Error ^ ^©-Errof ^ "J^-Error 
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0 4 6 8 10 12 

Number of Iterations 


niph J.l : ^ J in I>c«Blurring by De La Vall6c-Poussin integrals (Ti weighted image) 



Number of Iterations 


:;r«ph 32 : % -Mrmr in I)c- Blurring by De La Vall^c-Poussin integrals (Tz weighted image) 



Gn.pl. 3 J: » -En™ in D.-B.umng by D. U V=ll&-P«n i„.«g.bls (Dn.si.y weighted image) 
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A Iwk at Table 3.1 reveals that the %-Li-errors suddenly drop from 13.51 in 
the case of blurred image to its about 81% value 11.06 right in the case of the first 
iterate. The subsequent decrease of the errors is rather slov.’ but monotone. The final 
error H.57 in the 10th iteiate i.s about 63% of the initial error. The corresponding 
percentage drops in L.'- errors starting at the value 17.75, are from 79% to 55%, and 
for the L. eiror. starting at 46.26, the decrease is from 75% to 50%. 

In Table 3.2 eoriesponding to the T 2 -weighted image, the respective figures in 
the Li'Case are: 12.00. 88%, and 73*^; in the L;-case; 16.66, 88%, and 69%; and for 
the L..case: 65 HO. 83% , and 54%. 

In Table 3.3. coticsponding to the proton density weighted image, these are - 
in the I, i -case: 8.60. 81% , and 62%; in the L^-case; 11.39, 79%, and 56%; and in the 
L.-case: 42.3 \ 85%. and 71%. 


'fhese values arc suinuKiri/ed in the following Table 3.4. On the basis of the 
initial crrois it gets reiterated that the density weighted images are the smoothest, 
followed by the T: weighted and the Ti weighted images, which are approximately of 
the same srnoruhness, 'ITie exception occurs in the Ta weighted Loo-case where the 
initial error falls out of the line. Thm, as remarked earlier, is attributable to the 
generally fluctuating nature of the discretely computed Uo-errors. 



in itj,* I 1 fttoi fit’ni 
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tiutial 
% -bnof 

1 J.'hf 

heratf 

i,.iht 

Iterate 

tniiMl 
% -Error 

l-ir.st 

Iterate 


Initial 

%-Error 

First 

Iterate 
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Iterate 

"t, 

1T‘ 

! 

P i 


81% 

bn 

17.75 


55% 

46.26 

75% 

50% 

12,00 


Wi 

Upgl 

88% 

69% 

64.80 

83% 

54%. 

H.69 

''5i% 

62% 


SB 

56% 

42.33 

85% 

imiiiiii 


T.bk .1.4 • of Dc-Blumog of T„ T,. & Donsitj woightod images by Vallfe-PoKsmimegraJs 

'Hie intcresliny point ui be observed here is that even though the numeri 
em.rs for iterates do not reduce very subslamially, the visual effect of their de- 
blurring as well their visual nearness to the onginal is unmistakable. 





























3.4.3 SIGNED ERRORS IN THE ITERATIVE DE-BLURRINGS BY V'"' 
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Ihc iiitcirnwtion t'n the point-wi.se behavior of the convergence of the iterative 
de-bhii rings by Vj"' is presented in the following strip-graphs (Graphs 3.4 (a-b)). As 
in the previous chapters, the signed errors / - Vj"”/ in the reconstructions Vj™’ / of 
/ arc shown in the form of the strip-graphs for the first Ti-weighted cross-section. 
The errors arc displayed for the same horizontal 128”’ mid-row and for all the 256- 
pi \cls ini it. As in the case of strip-graphs in Chapters 1-2, the unit for these graphs is 
one gray Icud when the original is assumed to be a 256-gray level image. 

nic signed enors for the blurred image / ~T„/ s / -Vf^f are shown in the 
first strip graph. Tliesc errors for the de-blurred cases of the first ten iterates 
/ " m S 11) (Kcur in the next ten strip-graphs. The last strip-graph is a 

superposition of all the .signed errors in the ten iterates, as well as in the blurred 
image. 

As piiintcd out earlier in the case of the Fej6r sums in Chapter 1, here, too, 
some of the points where the values undergo an abrupt change have not got depicted 
in the strip graphs. This is indeed so for the cases m = 1, 2, and 3, where the range 
exceeds the limit 30 depicted in the graphs. 

These strip-graph.s show a similar, more or less monotone, behavior of the 
point-wise errots with increa.sing m. The errors, along the blackened raid line, may be 
observed tt> be relatively small for the pixels where the cross-section is relatively 
smooth and large where the cross section is less smooth. The error patterns for 
different m’s are similar where the magnitudes of errors keep decreasing as m 
increa.ses. A surge in the error activity, analogous to that seen with the earlier cases of 
Fejtfr and the generalized Jackson operators, for pixels near the boundaries of 
anauraical .structure is observable here too. Similar to what has been the case with the 
percentage errors presented in the tables before, the signed errors for the De la Vall6e- 
Pou.s.sin blurring and de-blurrings closely mimic the case of the generalized Jackson 
operators with p = 3, and n = 45. 
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CHAPTER 4: APPLICATIONS TO MRI DIAGNOSIS 


In tins final chapter we ronnd-up the work of the previous chapters by discussing 
sonic aspects of itie applicabilitv of the de-blurring scheme of Rathore, with a special 
emphasis ini his chtucc of hejcr sums as an easy to implement, general purpose blurring 
mmlcl for the routine MRI work. The simulations presented in the previous chapters 
clearly lienionstrateti the apcilicahiliiy of the scheme to improved visual diagnosis on 
actual MRI cross sectiotis. The de-blurrcd image.s, even with only moderately close other 
choices of blurring ttHHlels. arc relatively sharp, possess improved tissue contrast and 
boundary delineation, anil do not sccni to overemphasize the noise. 

In the following, we nunierically analyze and experiment with some of the 
questions and issues jxn taming to: (i) a practical choice of the system blurring operator, 

(ii) the tvliavior of the dc-hlutTing .scheme with respect to varying sizes of machine data, 

(iii) the mechanism of dc bhiiiing, and (iv) an application to improved numerical tissue 
parameter estimation, 

A miKlei of the system blurring i.s indeed the most crucial and important input to 
the scheme of dc blurring 'Oie question is whether the scheme is robust with respect to 
the choice of the blurring model, or does it too critically depend on it. The estimates of 
Rathore in C'haptcr 0 asserted that if j|dij| = ||I->4[1 < 1, 

AdVf)-f-^A-\V-A)f. 

and P«Ali<({I-Mr)/(l-IHI)]|l^l- 

It follows that unless the system blurring operator U is too far away from the 
identity operator I, any reasonable choice of A lying, so to say, in between the 17 and the I 
and more towards V in the sense of |117-I{| < ||A|1 1|I-A||, are expected to result in a 
reasonable dc-biurring. Moreover, since A« acts on the diffused noise A, unless both m 
and }|4|j arc large the noise is not expected to reassert itself. 
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4.1 A COMPARISON OF THE FILTER COEFFICIENTS 



<;rn|»h 4.1 . Filter (.’tHtfficients of Fej^r, Jackson (p = 2, 3, 4), and Vallde-Poussin Operators 


'I'he blurring operations considered in the Chapters 1-4 have been 2D- 
convolutiims with the ID-kernels of the type: 

= + coskt, 

where the coefricient.s of the kernels associated with the Fej6r operators, the generalized 
Jackson operators and the Vallde-Poussin operators are given as: 


Fejer operatttrs: pj"’ ~\-kln, 

Jackstui operatt'rs: p]"* = / Po where p^"”^ = ^(“1) 

(n!)' 


j-O 


. A 


np + p — n j-k-l 
2p-l 


,& 


VaIl&-Poussin operators: Pt 




(n~k)lin + k)l 
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The graph above (Graph 4 . 1 ) shows the filter coefficients, appropriately 
renumhcied and syninietrically folded, for the operators that have been used in Chapters 
I -3 for the simulaticms tliere. It makes clear that a gradation of these operator, according 
to an increasing order of (he associated bluiring, is given by: (t„, L„2, L„3, T„, and 1^4. 
The closeness of the graphs of filter values of Uy, and is clearly responsible for the 
near equal K'havior of these operators in all the simulations. 

On the basis of the above graph the Fejdr sums cr„ are expected to produce the 
least hluuing. and so the associated approximation errors ought to be the smallest for 
them. Amongst the Jackson opt‘rators the errors are expected to be the least for Uy, 
fidU'wed by i.r. u and I.„4. The Vallde-Poussin operators V„ should fall in between L„.3, 
and / ‘Hic following error tables, w’hich present the %-Li, %-L2, and %-Uo errors re- 
assorted from the tables of the previous chapters, readily illustrate this point: 
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4.2 AN EXPERIMENT WITH THE CROSS-FILTERING 


Indeed one of the principal advantages of the iterative de-blurring is its interactive 
monitoring feature of continuing the iterations only till the noise does not start reasserting 
itself. Practical experience with the iterates tells that each of these has something useful 
to tell about the cross-section studied. Quite often the individual iterates suggestively 
point out to some typical feature of the cross-section which subsequently gets accepted or 
rejected by a closer look at the original and all the available iterates jointly. The present 
section consists of simulations done in order to obtain an insight into this applicability of 
the de-blurring in the actual practice where the knowledge of the system blurring operator 
could normally be obtained only indirectly. Thus, if nothing better is available, one could 
start with a reasonable set of empirical filter functions and examine their applicability by 
obtaining a few trial reconstructions. Reconstructions based on well-calibrated phantoms 
may also be useful in generating approximations of the system blurring operator. 

In the simulations below the Ti, T 2 , and p-weighted cross-sections in Fig. 0.4 are 
blurred by each of the three operators ct„, Lnj, ip = 3), and Vn- As the generalized Jackson 
and the De la Vallfe-Poussin operators induce similar but a greater blurring than the Fej6r 
sums, the numerical experiment below provides insight into the nature of de-blurring 
when the assumed blurring model results in a slightly more (or less) smoothing compared 
to that of the actual system blurring operator. Thus by modeling the system blurring 
operator after the above mentioned operators On, Lnj, ip = 3), and V„, each blurred image 
is subsequently iteratively de-blurred by the associated de-blurring operators. The 
resulting images may be readily cross analyzed. The experimental findings point out that 
the de-blurrings associated with the Fejer operators (7„^m result in sharpened images of 
least error out of all of the operators considered here. This methodology, with appropriate 
models of the blurring, could be used empirically for estimating system blurring on an 
MRI machine. The results also indicate that neighboring models perform more or less 
equally well and thereby point out to the robustness of the de-bluning technique. The 
images and tables presented in Fig. 4.1-4.3 and Tables .4-4.6 correspond to the third (m 
= 4) iterates, which in actual practice have been found to be quite satisfactory. 



110 





Fejer- Jackson (p=3) 


Fejer- De La Vallee-Poussin 


JatLson -Fcjci 


Jackson (p=3)-- Jackson (p=3) 


Jackson (p=3) -De La Vallee- 



Ik La Fej<fr 


Vall6e-Poussm-Vall6e-Poussm 


Fig. 4 1: Blurring and De-Blurring by Different De-Blurring Operators on a Ti-Weighted Image 


1 

Fej6r 

Jackson (p=3) 

Vall6e-Poussm 


5 501431 

6 843048 

7.142781 


10.955830 

9.755668 

9.755668 

1 V'allec-Poussin 

11.143796 

9.882319 

9.836020 


Li Error 



Fejdr 

Jackson (p=3) 

Vall^e-Poussin 


5.805688 

7.658533 

8.092607 

Jackson(p=3) 

13.614547 

11.791772 

11.791772 

ValMe-Poussin 

13.903517 

12.052172 

11.938289 


La Error 


1 

Fejdr 

Jackson (p=3) 

Vall6e-Poussm 


14.031715 

16.998110 

18.222013 

Jackson(p=3) 

32.352612 

27.978689 

27.978689 

ValWe-Poussin 

33.391159 

28.647913 

28.262102 


L JError 


Table 4. 4: Blurring and De-Blurring by Different De-Blurring Operators on a Ti-Weighted Image 
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De Ijl Vallfe-Poussin- Fej<5r De La ValMe-Poussin-Jackson (p=3) Vallee-Poussm-Vall6e-Poussm 

Fig. 4. 2: Blurring and De-Blurring by Different De-Blurrmg Operators on a Tz-Weight Image 



Fejdr 

Jackson (p=3) 

Vall6e-Poussin 

Fejdr 

5.909182 

6.485254 

6.668396 

Jackson(p=3) 

10.502147 

9.730365 

9.717400 

Vallte-Poussin 

10.612084 

9.801308 

9.776073 


Li-Error 



Fei6r 

Jackson (p=3) 

Vall6e-Poussm 

Fejdr 

7.100663 

7.104949 

7.369535 

Jackson(p=3) 

14.437531 

13.189909 

13.139686 

Vallte-Poussin 

14.586150 

13.302277 

13.239385 


Lz-Error 




Jackson (p=3) 

Vallle-Poussin 

Fejdr 

19.793056 

15.321804 

16.175762 

Jackson(p=3) 

51.579433 

43.934147 

43 159855 

ValUe-Poussm 

52.697315 

45.431942 

44 509380 


L ..-Error 


Table. 4.5 : Blurring and De-Blurring by Different De-Blurring Operators on a Tz- Weighted image 







112 





Dc Valltfe-Poussin- Fejt*r 
Fig. 4.3: Blurring and De-Blurring by 


De La Vall^e-Poussin-Jackson (p=3) Vall6e-Poussin-Vall6e-Poussin 
Different De-Blurring Operators on a Density- Weighted Image 



Fejdr 

Jackson (p=3) 

Vall6e-Poussin 


3.477801 

4.455467 

4.660405 

Jack.sim(p=3) 

7.061089 

6.221765 

6 206601 

Vall^e-Pou.ssin 

7.187753 

6.307024 

6 274099 


Li Error 



Fei6r 

Jackson (p=3) 

Vallte-Poussin 

Fei6r 

3.828713 

5.043094 

5.337027 

Iack.son(p=3) 

8.760753 

7.623755 

7.574858 

Vallfc-Foussin 

8.950841 

7.778380 

7.707120 



Fei6r 

Jackson (p=3) 

Vall^e-Poussin 

Fei6r 

15.780841 

14.786440 

15.769584 

.... * 

Jackson(p=3) 

34.737602 

33.660534 

33.788979 

Vallfc-Poussin 

35.423737 

33.459160 

33.582581 


Lo, Error 


Table 4.6 : Blurring and De-Blurring by Different De-Blurring Operators on a Density-Weighted Image. 
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4.3 NUMERICAL RESULTS ON DE-BLURRINGS AT 
DIFFERENT RESOLUTIONS 

The theoretical results in the previous chapters were concerned with the behavior 
of the de-blurring operators as the parameter n, related with the data size there, 
approaches infinity. Even though the results predict improvements in the reconstructions 
with increasing n, in the actual practice we are constrained to work with the moderate 
values of n such as 256. In this connection, the theoretical results in conjunction with 
realistic numerical .simulations, nevertheless, help us in making projections about the 
future practice when the larger sizes may become possible due to improvements in the 
technology. Also, realistic numerical simulations provide an added incentive for the 
needed technological developments if the results at the larger sizes are expected to be an 
order of magnitude better, as compared with the present state of the art. 

The next experiment in this chapter is intended in the same direction. In it we 
examine the efficacy of the de-blurring on different sizes of cross-sectional data, or 
equivalently on the MRI images of different resolutions such as 64x64, 128x128, 
256x256, 512x512, etc. Much of the present clinical practice is based on 256x256 
machine data. With the discovery of newer and faster pulse sequences and improved 
machination, it is expected that soon in the future routine MRI diagnosis may switch over 
to the 512x512 case. Images at the 512x512 resolution would be capable of revealing 
finer details than is possible with the 256x256 images. Indeed it is possible right now to 
acquire the 512x512 data. However for it to match the noise characteristics of the 
256x256 images the required imaging times are prohibitive at the moment. 

To carry out a realistic numerically-simulated error analysis, we used a Tr 
weighted 5 1 2x5 1 2 spin echo image data at TE = 15 ms and TO = 400 ms, to generate the 
“Original” cross sections of sizes 64x64, 128x128, 256x256, and 512x512, from which 
the respective bluired images were obtained by using the Fej6r filter coefficients for the 
size 512x512. This is keeping in mind that the basic machine out-put is this data for the 
image and not the image itself. As an implication, the blurring would be maximal for the 
image of the size 512x512 and minimal for that of 64x64. The de-blurrings were done by 
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the corresponding Fej6r filters. Two kinds of errors were calculated: (i) on the 512x512 
grid of the original, in which the lower size images were extended to the size 512x512 by 
using zero padding in the frequency domain (the extended images thus representing the 
values of the same trigonometric polynomial on a larger grid of points), and (ii) the 
indi\'idual errors in reconstructing the images themselves of the respective sizes 64x64, 
128x128, 256x256, and 512x512. Thus the percentage errors in (i) are with respect to the 
“Originaf’ of size 512x512, and in (ii) these are with respect to the “Originals” of the 
sizes 64x64, 128x128, 256x256, and 512x512. One expects large errors in (i), and much 
smaller eiTors in (ii), ,since the filtering is progressively less for the smaller image sizes. 
Hence, in the case (i) the errors are not supposed to decrease any significantly, they being 
predominantly due to the basic errors in the approximation by the partial sums of the 
Fourier series associated with the image sizes 64x64, 128x128, 256x256, and 512x512. 
However, these errors for the sizes 64x64, 128x128, and 256x256 in case (i) do give us a 
rough idea of the actual errors incurred in MRI reconstructions with the respective sizes 
64x64, 128x128, and 256x256. The errors in both (i) and (ii) for the size 512x512 are 
indeed to be the same and they reflect the reduction of the errors in the successive de- 
bluning iterates, as they do for the remaining sizes in (ii). 

The related images are shown in Figs. 4.4-4.7. The numerical errors in the case-(i) 
are presented in the Tables 4.7-4.9 and the same are also shown in the Graphs 4.2-4.4. 
The same in the ca.se-(ii) are given in Tables 4.10-4.12 and Graphs 4.5-4.7. As seen from 
these, the nature of the errors is, indeed, as described in the previous paragraph. The 
agreement of the order of the errors in Tables 4.10-4.12 for the size 512x512, and in 
Table 1.1 (for the Ti-weighted cross-section there in Chapter 1) seem to confirm the 
appropriateness of the model of the filters chosen. On the other hand, the relatively too 
small errors in Tables 4.10-4.12 for the sizes 64x64, 128x128, and 256x256 indicate a 
better performance of the de -blurring iterates if ||A|1 (as described in Chapter 1) is 
relatively small. 

We also include a repeat of experiment of Fig. 1.1 in Fig. 4.8, where the 512x512 
image is regarded as a machine output (with 512x512 Fej6r blurring model)and the de- 
blurring iterates and the result of a direct filter inversion are shown. The improvements 
are as before and a similarity in the results of direct inversions in the two cases is striking. 
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Iteration I Iteration 4 



Iteration 7 Iteration 10 


Fig. 4.4 : De-Blurring by Fej6r sums (64x64 Tr Weighted Image) 
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Iteration 7 


Iteration 10 


Fig. 4.5 : De-Bluning by Fej^r sums (128x128 Tj-Weighted Image) 
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Iteration 1 Iteration 4 



Fig. 4.6 : De-Blurring by Fej6r sums (256x256 Ti-Weighted Image) 
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Blurred 


Iteration I 


Iteration 4 


Iteration 10 


Itemtion 7 


Fig- 4J : De-Blurring by Fej6r sums (512x512 TrWeighted Image) 
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256 

512 

12.760939 

1 1 232959 

11.327273 

8 407426 

11.159102 

7.322424 

11 127623 

6.648437 

11 122073 

6 161029 

11.122254 

5.782414 

11.123207 

5.475666 

11 124063 

5 219060 

11 124699 

4.999608 

11.125145 

4 808920 

11.125442 

4 640935 


14 968604 


14.984726 
14.990157 
14 991931 
74.992516 
14.992712 
14.992780 
14.992803 
14.99280" 


I able 4.7 ■ “a I .j I'nors in recon.structing original 512x512cross-section by image data of different 
M/cs (bhitTing & de-blurring.s associated with Fej6r sums for n = 256). 
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26.063387 


1'ablc 4.8 ; 1,. -Error lirror.s in reconstructing original 512x5 12cross-section by image data of different 
sizes (blurring & de-blurrings associated with Fej6r sums for n = 256). 
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69.105209 
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68.0.34599 



128 

256 

512 

16 756004 

12 835905 

1 1.379604 

15 765809 

10.638448 

7.717406 

15.709594 

10.357162 

6 543892 

15 703228 

10 293384 

5 875368 

15.702043 

10 274142 

5.410994 

15 701744 

10 267305 

5 058268 

15.701655 

10.264592 

4.776148 

15 701626 

10 263425 

4 542590 

15.701616 

10.262887 

4 344323 

15.701613 

10.262620 

4.172835 

15.701612 

10.262481 

4.022374 



128 

256 

512 

53.528534 

41.391506 

37.572300 

48.251698 

33.933582 

26.652369 

46.826370 

32.256283 

21.929525 

46.472794 

31.923101 

18.960024 

46.471657 

31.895712 

16.670185 

46.468437 

31.924402 

14.764875 

46.466652 

31.952774 

13.141476 

46.465881 

31.971827 

11.748270 

46.465652 

31.983000 

10.899012 

46.465771 

31.989004 

10.190278 

46.466305 

31.991884 

9 577790 


i 68.0.34599 | 46.466305 | 31.99i»b4 | 

Table 4.9 : L. -Error Errors in reconstructing original 512 x 5 I 2 cross-section by image 
sizes (blurrint! & de-blurrines associated with Fe.6r sums for n = 256) 
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Number of Iterations 

(Jraph 4.2 : L- l'iu>t I'rtors in reconstructing original 5 12x51 2cross-section by image data of 
thf fcH’iit M/cs ihluniiu’ k tk-Wurrincs associated with Feidr sums for n = 2561. 
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Number of Iterations 

Jraph 43 : Lj-F,rrur linors m reconstructing original 512x512cross-section by image data of different 
sizes fWnniiic k ik blurrines associated with Fekr sums for n = 256). 



Number of Iterations 


Graph 4.4 : L. -Error liintrs in reconstwcling original 512 x 5 l 2 cross-section by image data of 
different sizes f hlumne & de-blumnes associated with Fei6r sums for n = 256). 
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2 418032 


0,324249 


0 051506 


0.009190 


0 002316 


0.001470 


0.001450 


0.001500 


0.001622 


0.001912 


0.002547 


128 


3.818776 


0 860875 


0242072 


0.075348 


0025235 


0.009099 


0.003753 


0,002055 


0.001666 


0.001839 


0002461 


256 


5 611177 


2 270876 


1 191611 


0 689982 


0420825 


0.265632 


0.171963 


0 113500 


0076070 


0 051608 


0 035410 


512 


11.232959 


8.407426 


7.322424 


6.648437 


6.161029 


5.782414 


5.475666 


5 219060 


4.999608 


4.808920 


4.640935 


Table 4.10 : %-U Errors m reconstructing onginal cross-sections of different sizes (blumng & de- 
blurrings associated with Fejdr sums for n = 256). 
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128 

256 

512 

3.420601 

5.804968 

7.707932 

11 379604 

0.432999 

1.271163 

2.765119 

7.717406 

0.067814 

0.346227 

1 351323 

6.543892 

' 0.012200 

0.105967 

0 754756 

5.875368 

0.003433 

0.035028 

0.451532 

5.410994 

0.002680 

0.012368 

0.281519 

5 058268 

0.002748 

0.004987 

0.180475 

4.776148 

0.002853 

0.002930 

0.118067 

4.542590 

0.002965 

0.002584 

0.078458 

4.344323 

0.003140 

0.002669 

0.052805 

4.172835 

0.003559 

0.003091 

0.035959 

^ 4.022374 



Table 4.11: %-L 2 Errors in reconstructing onginal cross-sections of different sizes (blurring & de- 
blurrings associated with Fejbr sums for n = 256). 


512 
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18.960024 
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128 

256 

7.346234 

13.665717 

24.851616 

0.979912 

3.829877 

10.249213 

0.151225 

1.180119 

5.414123 

0.024498 

0.382845 

3.110082 

0.008779 

0.128460 

1.884815 

0.009864 

0.043945 

1.182202 

0.010758 

0.015377 

0.758503 

0.011491 

0.009114 

0.494281 

0.012109 

0.009488 

0.325697 

0.012641 

0.009688 

0.216301 


0.013116 


0.009726 


0.144287 


I ■ 

Table 4.12 : %-L« Errors in reconstructing original cross-sections of different 
blurrings associated with Fej6r sums for n = 256). 
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Number of Iterations 

Graph 4.5 : %-L{ Errors in reconstructing original cross-sections of different sizes (blumng & de-blumng: 

associated with Fejdr sums for n = 256). 



N umber of Iterations 

Graph 4.6 : %-L: Errors in reconstructing original cross-sections of different sizes (blurring & de-blun 

associated with Fejer sums for n = 256). 
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Graph 4.7 : %-lL. Errors in reconstructing original cross-sections of different sizes (blurring & de-blunings 

associated with Fej6r sums for n = 256). 
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4.4 AN APPLICATION TO IMPROVED TISSUE 
PARAMETERS ESTIMATION 


The next numerical experiment in this chapter explores the use of iterative de- 
blurring in obtaining a relatively more accurate set of Ji, T 2 and Proton density (p) 
weighted images from the corresponding blurred images, a use of which in the 
calculation of the corresponding T\, Tj and Proton density maps may reasonably be 
expected to give better results. 

A calculation of Ti, 72 and Proton density (p) values indeed involves the solutions 
of the Bloch equations for the various imaging protocols associated with different pulse 
sequences. For the SE images, e.g., one uses the equation 

where the function / denotes the image intensity and the rest magnetization M is 
proportional to the proton density p. Thus, given a set of independent cross-sectional data 
using three appropriate sets of EE’s and r/?’s, the tissue parameters T\, T 2 and proton 
density (p) could be obtained by solving the related system of equations for each pixel, 
where in principle any convergent numerical solution procedure is supposed to give the 
same values (we have used an in house method developed by S.B. Rao et al.(2000)). 

We utilized the Ti, T 2 and Proton density weighted images of Fig.0.4 and tissue 
parameters were calculated from these as described before. These were taken as the 
reference Ti, T 2 and proton density maps of the cross-section, using which, Ti, T 2 and 
density weighted spin echo images at the same pulse sequence parameters of the original 
images (i.e., TE/TR = 15/300, 20/20(X), 80/2000 ms) were synthesized. These image sets 
were blurred and iteratively de-blurred by the Fej6r sums, and the tissue parameter values 
were computed from each of the reconstructed image sets. . The resulting Ti, T 2 and p 
maps are shown in Figs. 4.9-4.11. From these the associated percentage Li-errors (with 
respect to the original reference maps) have been computed and are presented in Tables 
4.13 and depicted in the Graphs 4.8. The resulting improvements in the tissue parameter 
calculations with the successive de-blurrings are obvious from these. 
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Ongmal Blurred 



Iteration 7 Iteration 10 


1 


Fig. 4.9 : Ti maps calculated from synthesized images blurred and iteratively De-Blurred by Fejer sun 
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Iteration 4 



Ilrfiitiiitt 1 


Iteration 10 


Fl«. 4.19 ; T, imp, mklitaol from s) nlhcw.^a imaeo.r Wurred and ileralively De-Blunnd by Fejdr sums 
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Iteration 7 


Iteration 10 


Fig. 4.11 ; p maps calculated f»m synthesised images blurmd and itendvely De-Blu™d by Fejdr 
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1’abk* 4.13 : ';?-Li I-iiTors in Tl, T 2 and p maps calculated from synthesized images 
bluned and iteratively De-Blurred by Fejdr sums 


— . — T1 



Graph 4.8 : %-U Krror.s in T,, T, and p maps calculated from synthesized images blurred and 

iteratively De-Bluired by Fejdr sums. 
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Finally, we conclude the work by appending: (i) the graphical output for the discrete 
point spread functions associated with the blurring and de-blurring operators considered 
in the previous chapters, and (ii) the error maps (considerably brightened to show the 
pattern of the errors) associated with the de-blurred reconstructions of the Ti-weighted 
image (Fig. 0.4) for all the operators considered previously. These, so to say, provide the 
information on “de-blurring in action” and pictorially illustrate the working of the de- 
blurring scheme. 

The first graph below shows the psf s of the Fejdr sum and the related de-bluning 
operators. We note that whereas the central output of the Fejdr sum lies between 0.4- 0.6, 
the out put of the 10* de-blurring operator is near 0.9. Moreover, the psf of the Fejdr sum 
contributes only positively (leading to blurring), the negative contributions of the de- 
blurring operators are responsible for the de-blurring by them. The behavior of the 1“ to 
the 9* de-blurring operators is in between that of the Fejdr sum and the 10* de-blurring 
operator. The psf s of the generalized Jackson and the De la Vall6e-Poussin operators 
follow next. 



Graph 4.9 . PSF’s of De-Blurs (by Fej^r sums) 
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Graph 4.10 • PSF’s of the De-blurs ( Generalized Jackson operator (p=2)) 
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Graph 4.11 : PSF’s of the De-blurs ( Generalized Jackson operator (p=3)) 
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Graph 4.12 : PSF’s of the De-blurs ( Generalized Jackson operator (p=4)) 
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Graph 4.13 : PSFs of the De-blurs ( La ValI6e-Poussin) 


We note, from the Graph 4.10-12 of the psf s related with the generalized Jackson 
operators, that in the respective cases of p = 2, 3, and 4, the central outputs of the basic 
operators and the last of the de-blurring operators are respectively between 0.3-0.4, and 
near 0.7 for p = 2; near 0.3, and near 0.6 for p = 3; between 0..2-0.3 and between 0.5-0.6 
for p = 4; and the same for the De la Vallde-Poussin are near 0.3, and near 0.6 ). The 
positive and negative output character of all these cases is similar to the case of the Fejer 
sum. The central outputs in all the cases increase monotonically in value. These graphs at 
once corroborate the earlier observed approximation and the blurring and the de-blurring 
properties. 

Thus, the peaks near .5, .4, .3, below .3 and near .3 of the basic operators grade 
their blurring properties, whereas those near .9, .7, .6, below .6, and near .6 of the last 
iterates grade the corresponding order of approximation properties of the associated de- 
blurring operators. Finally the sign changes observed for all the de-blurring operators 
assert their de-blurring attribute. 
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Fig. 4.12; Magnitudes of Difference Images in De-Blumng by Fejer sum 
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Original-Iteration 7 



Original-Iteratian 8 Ongmal-Iteration 9 Onginal-Iteration 10 


Fig. 4.13 : Magnitudes of Difference Images in De-Blurring by Generalized Jackson Operator (p=2) 





135 



Onginal-Iteration 2 Onginal-Iteration 3 Ongmal-Iteration 4 



Original-Iteration 5 Original-Iteration 6 Onginal-Iteration 7 



Onginal-Iteration 8 Original-Iteration 9 Onginal-Iteration 10 


Fig. 4.15 : Magnitudes of difference Images in De-Blurring by Generalized Jackson Operator (p=4). 
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Ongmal-Iteration 2 Onginal-Iteration 3 Ongmal-Iteration 4 



Ongmal-Iteration 5 Onginal-Iteration 6 Ongmal-Iteration 7 



Ongmal-Iteration 8 


Ongmal-Iteration 9 


Ongmal-Iteration 10 


Fig. 4.16 : Magnitude of Difference Images m De-Blurring by De La Vall6e-Poussm Operator 
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The intensified images of the errors presented in Figs. 4.11-4.16 reiterate the 
qualitative behavior suggested by the graphs of the discrete psf s, discussed before. From 
these we also observe the progressive increase in the error of approximation by the de- 
blurring iterates for the Fejer sums (n =128), the generalized Jackson operators (with p = 
2, n = 65; p = 3, n = 45; and p = 4, n = 35), respectively, and the near similarity of the 
error images for the case of the generalized Jackson operators (with p = 3, n = 45) and the 
De la Vallde-Poussin operators (n = 128). The interesting common feature, however, has 
been the visually unmistakable de-blurring property (in all the reconstructions presented 
in Chapters 1-3, including in the case of generalized Jackson operators for p = 4, in spite 
of their poor approximation property) observed for all the operators subjected to de- 
blurring scheme of Rathore. 

In conclusion, the approximation and the de-blurring properties are two distinct 
aspects of the same de-blurring operators 

”C/^-’,m = l,2,3,4,... 

which the theoretical and the simulation results of the thesis have tried to clarify and 
illustrate with the operators a„, Ln^, and V„, considered for the purposes of application in 
the field of MRI diagnosis. 
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